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SUMMARY 

This document presents methods for developing the ordinary differential equations 
(o.d.e.) of motion in canonical form equivalent to the forms of Delaunay and Poincare. It 
also presents modifications to these forms so that three variables, which are constants of 
motion, result while the forms remain canonical. 

The equations of motion are for a vehicle propelled by constant thrust magnitude with a 
constant mass flow rate. The vehicle is moving in a central force field. The trajectories 
are optimum in the sense of classical calculus of variations in a neighborhood definable by 
the boundary conditions of the specific problem. Specific problems are not discussed in this 
document. 

The value of the document lies in two major areas: 

1. The possible economics in numerical calculations which may result from using 
these ordinary differential equations, and 

2. The application of the general perturbation theory of classical celestial mech- 
anics to approximate solutions of these ordinary differential equations. 

This document has been written to record the results of the investigation and was not 
meant to be a tutorial treatment of the subject. For such treatment, references 1 through 4 
are recommended by the author. 
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2. Goldstein, H. : Classical Mechanics . Addis on -Wesley Publishing Co., Inc., Cambridge, 
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3. Ford, L.R. : Differential Equations. McGraw-Hill Book Co., Inc., N.Y., 1933. 

4. Smart, W.M. : Celestial Mechanics . Longmans, Green, and Co., Ltd., London, 1953. 



I. INTRODUCTION 


This report addresses itself to the development of the ordinary differential equations 
(o.d.e.) governing the motion of a vehicle propelled by constant thrust magnitude with con- 
stant mass flow rate in a central force field. The trajectories are optimum in the sense of 
classical calculus of variations. The variables used to define the trajectory are canonical 
and based on the solution of the trajectory when the thrust magnitude is zero. 

The Hamiltonian for the full problem is given by 


where 


H = K + 


w 


H 3 T 3 2 2 

r r cos e r 


2 

- x 2 ^2" Sec 9 tan 9 + p i u 

z r 


+ p 2^2 + p 3^ SeC * 9 “ H A(X) 

r r 


o- = -m 
u = r 


v/r = $ 
2 2 

(w/r ) sec e = 4> 


A(X) = 




,2 2 ,2 2 2 „ 2 

' ■ - ' + r cos e x 3 


and X j, \ 2 > V 3’ p i> p 2’ p 3 ’ and \ are the Lagrange multipliers. 


( 1 ) 


We consider here 

H = H q + H x (2) 

F 

where = — A ( x ), and we present various solutions of H q (base problem). 

This report proposes only to present methods of deriving the various transformations 
and differential equations. Only enough detail is given so that the results may be reproduced. 
This paper is concerned only with methods. The basic theory may be found in the references. 

The approach was based on solving the problem for the Delaunay variables, and then 
transforming the Delaunay variables to other variables. This was accomplished with ease. 
However, the inverse sets of transformation^ created problems in defining the new variables 
in the initial coordinates of (r, 9 , <j>, u, v, w, X, p, X^). These were algebraic problems 
and were solvable. However, this does not assure a correct answer. Therefore, a more 
direct approach was sought. Noting that these transformations form a group, one is assured 
of a direct transformation from the initial coordinate systems to the final desired coordinate 
system. The following diagram illustrates the problem: 
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In the diagram, the dashed lines represent the first approach and the solid ones the 
approach used in this document. 

The problem given by Eq. (1) has four constants: the Hamiltonian and three others. If 
these last three constants are made to be coordinates in the transformed sense, then there 
are three less integrals to evaluate in numerical integration or in perturbation procedures. 
Therefore, in two of the sets of variables these are made coordinates. 

The author desired one transfer function in which the constants were to be modified for 
all four cases. This was not accomplished. However, two transfer functions did do the job 
except that for the modified Poincarg variables added transformations are needed. The 
Delaunay and Poincare variables are on state only. One transfer function, based primarily 
on state, serves for this type of transformation. The cases with the three constants were 
dependent for some variables on both state and Lagrange multipliers. One transfer function, 
with additional transformations for the modified Poincare variables, serves for this type of 
transfer mation . 

The approach used shortens the path to solutions of the problem but does not eliminate 
much detailed algebra and differentiation. Therefore, support for this type of work was 
obtained by use of FORMAC. * It was decided that FORMAC could be used favorably for the 
following steps in the development: 

1. Obtaining the partial derivatives of the generating function, with respect to the 
coordinates Q and substituting them into the Hamilton -Jacobi partial differential 
equation to confirm the generating function; 

2. Obtaining the first partial derivatives of the generating function with respect 
to the transformed momenta [K (a)] for the transformed coordinate L or (3; 

3. Obtaining the second partial derivatives of the generating function with respect 
to the coordinates and momenta to display information for forming the ordinary 
differential equations of the full problem. 

The above work was supplemented by hand calculations at points so that the programs 
could remain general and manageable. 


*There were many detailed discussions and close liaison between the author of this document 
and Mr. D. Valenzuela of IBM. His report on the FORMAC work is in preparation. Fur- 
ther information may be obtained from Mr. Valenzuela, IBM, 1730 Cambridge Street, 
Cambridge, Massachusetts. 
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The results of this work have the following potential uses: 

1. Development of approximate closed-form solutions for trajectory problems by 
using the methods of celestial mechanics perturbation theory; 

2. Study of the characteristics of types of trajectories using the methods of celes- 
tial mechanics; 

3. Development of computer programs for numerical integration of the ordinary 
differential equations. 

Because of the canonical form this formulation should be easily expanded to include 
variations of the gravitational field. Also potentially large steps may be made in numerical 
integration because of the form of the variables. 

The information in sections n and m has been treated separately because of the designa- 
tion of the variables as momenta or coordinates. This designation is different for the 
Delaunay and Poincare variables and the modified variables. The organization of the text 
follows: 

1. The Hamiltonian for the base problem is repeated. 

2. A transformation to designate convenient momenta and coordinates is made 
for this set of problems. 

3. The Hamilton- Jacobi equation for this set of problems is presented. 

4. An indicated procedure for solution is presented and the general generating 
function for the transformation is presented. 

5. The coordinates are then presented in transformed form. 

6. A discussion of the ordinary differential equations is given. 

7. The perturbation function is presented. 

8. The logic for selecting the K.(ofj) for the particular coordinates is given (either 
Delaunay or Poincare) and the generating function for the transformation 
presented. 

9. Steps 5, 6, and 7 are then repeated. 

The procedures for obtaining the information presented are not unique. The choice of 
the procedures used was based on the following: 

1. The use of FORMAC to minimize errors; 

2. The ease in presentation; and 

3. The overall logic of the approach. 
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The reader is cautioned that the standard astronomical notations of « and p are used in 
both Poincare and Delaunay computations. Thus: 



D 

P 

K 

a 

a 

L 

P 

P 


It will be obvious to the reader from the text what a and P represent. The same holds true 
for the modified Delaunay and PoincarS variables. 

GLOSSARY 



r = radius vector to vehicle 

8 = angle out of reference plane ( e max = inclination) 

4> = angle in reference plane 

u = momentum conjugate to r (r) 

2 ' 

v = momentum conjugate to 8 (r 8) 

2 2 

w = momentum conjugate to <t> (r cos 8 $) 
t = time 

<r = negative of mass flow rate (- m) 
k = gravitational constant 
m - mass 

X^ = Lagrange multiplier connected to u 
X^ = Lagrange multiplier connected to v 
Xg = Lagrange multiplier connected to w 
X^ = Lagrange multiplier connected to m 
P i = Lagrange multiplier connected to r 
Pg = Lagrange multiplier connected to 8 
Pg = Lagrange multiplier connected to <f> 

The above terms are basic to this document. Other terms are included, but since the 
thesis of this work is based on a series of variables, no attempt will be made to define each 
and every one because of the confusion it would probably create. 
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n. DELAUNAY AND POINCARE VARIABLES 


This section presents a method of deriving a generating function in terms of parameters 
(K) which are to be regarded as functions of the new momenta of the problem. The base 
two -body Hamiltonian for the trajectory problem may be written as: 


H 


h 


2 9 

w sec * 0 
~ 


- “TT - X. 


w 


2 2 
r 


sec 0 tan 0 


+ pj u + p 2 — 2 " Pg ~ 2 ~ sec 0 - Xi 


where 


7 a 


( 3 ) 


-m = constant 


V/ 3 


= r 


= 0 


(sec 0) = 4> 
r 

H’ *2’ p l’ P2 ? p 3 ? and ^7 are the Lagrange multipliers. 

The first step will be to transform this equation by the following transfer function: 

S = PjXj + P 2 *2 + p3*3 + p 4 r + P5P2 + P 6P3 + P 7 m ( 4 ) 

The results of this tranformation give the relationships: 

as 


Q i = v 

Q 2 = K 2 

Q 3 = K 3 
Q 4 =r 


Q* = Pc 


Q 6 = p 3 

Qr, = m 


" u = P 1 aQj’ 

_ 9S 

_V = P 2 ” 9Q 2 

_ 9S 

" W = P 3 = 9Q; 

~ -D - 8S 

P 1 ~ P 4 “ 9Q 4 
-e = P = J§- 

5 8 Q 5 

* ~ p - 8S 

* P 6 9Q g 

X 7 ~ P 7 " 9Q " 


( 5 ) 



where the partial derivitives are listed for ease in writing the Hamilton -Jacobi equation 
which follows: 


0 = 


3S,q 
W + Q 1 


3S 


^ / 


as 


Q, 


2 as 

sec 


k 


^2 


as 

W 

3, 




sec 


2 3S 


tan 


3S 3S 


3S 

" 3Q 4 8Q 


- Qc 


as 

aQ, 


2 q: 


sec 


2 as 


- Q 


as 

6 bq: 


3Q, 





(6) 


We note that Q„, Q 3 , and t do not appear in the equations, so that, by separation of 
variables, we have:' 


as _ 

w ~ 


Ki, 



as 

*5 


-K„ 


(7) 


This new partial differential equation has the following differential equations of the 
characteristic strip: 


as 



2 as 

sec 


Q, 


tan 


as 


as _ as l 

ay 5 q2 


as 



t _2 2 as 

+ K, sec - 

3 ^5 k 

^ 


«4 


as 

^6 


as 

2 

sec 



as 


( 8 ) 


(9) 


( 10 ) 

( 11 ) 

( 12 ) 
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where the dot indicates the derivative with respect to a parameter t 

/8‘S d 3 S \ 

V dT 3Q^ ) 

These equations are simply the two -body equations which we integrate as follows: 
From Eqs. (8) and (9), we obtain an integral. We substitute this integral into Eq. (10). 
Next, we combine Eqs. (10) and (11) to obtain a second integral. The next integral is ob- 
tained from Eqs. (9) and (12). The last integral is obtained by use of the solution of Eqs. 
(8) and (9) in combination with Eqs. (9) and (11). 

The solutions are: 



= ^ - 4 


_ 

9Q 



K rj = cos 


-J2 _ i r\ K.sin ar . 

K 4 kQ 4 . -1 4 8Q 5 

9 4v /* 2 -Ki<C S 'V K 4- k ! 


(13) 

(14) 


(15) 


(16) 


fpr- = sin"* sin i sin (K 7 ~ f) (17) 

9Q 5 7 

where 

, ^~ kQ 4 

C0S 

sin i = tp 

K 4 

The terms u, v, f , and i used in the above equations are to be considered as symbols and 
not variables. 

Eqs. (13) through (17) substituted into the Hamilton -Jacobi equation (6) yield an 
expression for 
8 4 



On integration one obtains: 


S (Q 4 ) = (K x - K 2 ) 


uQ 4 k 

^ ' 27T cos 


.-1 k ' K 5 Q 4 


^/ ] 


k 2 - K 2 K, 
4 5 


- Qj u - Q 2 v 


( 18 ) 


+ Qg sin 


_! sin (K ? - f) 


sm l 


+ Q c 


^1 - sin 2 i sin 2 (K^ - f) 
This equation was combined with those derived before to obtain a generating function 


-Kg + sin~ * 


cos i sin 


M 


K o 

S = Kjt + Q 7 + (Kj - K 2 ) 


uQ 4 k _ik-K 5 Q 4 

K r ' ZW C0S 


K 5 y k 2 - KgK* _ 


(19) 


XT* 

. Olr _ 1 

-Kg Qg - Qj /-Kg + q - J + Q5 sin t sin i sin (K 7 - f)] 

4 Q 4 


-Q, 


sin i cos (K^ - f) 


- sin 2 i sin 2 (K^ - f) 
where S is considered as S(K,Q, t) = S and 


+ Q c 


-Kg + sin 


^ cos i sin (K^ - f) 


L - sin^i sin^ (K^ - f) 


K x = - H o , K 2 = a- \ 7> Kg = w, K 2 = V 2 + w 2 sec 2 e 


Kt 


K =-u 2 + ^-“ 4 
5 r r 


— j , Kg = <t> - sin 


Kg tan 6 


K 2 K 2 
K 4 • K 3 


(20) 


K 7 = cos 


-1 K-kr 


- sin 


K^ sin 9 


k: 


~ K 5 K 4 JA- “3 

The first partials of S with respect to K^ produce the conjugate variables 


9 S 

( i = 1, 2.... 7)o These are presented below and information to reproduce 

them will be given in a document to be published by IBM (see footnote on page 3). 
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Q.ii 

L i = t+ -jr- 

1 K 5 


k -1 k-K 5 Q 4 

COS — ; - — 


K 


3/ 2 


k 2 -K 2 K 
4 5 


Q 7 

L 2 "^ +<t - L l> 


L 3 - Q 2 


cos icos (K^ - f) 


— - Q P 


cos i sin (K^ - f) 


sin i 


1 - sin 2 i sin 2 (K^ - f) 3 ^ 2 sini Jl- sin 2 i sin 2 (K ? - f) K 4 


+ Qc 


sin (K^ - f) cos (K„ - f) 


-Q. 


K , 


2 2 
1 - sin i sin 


(K ? -f) 


L 7 ^2 


2 

sin i cos i sin (K^ - f) 


1 - sin 2 i sin 2 (K^ - f) 


3/2 


+ Q c 


sin icos (K^ - f) 

1 - sin 2 i sin 2 (K^ - f) 


+ Q c 


cos i 


1 - sin 2 i sin 2 (K^ - f) 


( 21 ) 


L 6 ” _Q 6 


L 5 - Q 1 ~2u + 


1 K 4 L 7 (K^-kQ 4 ) 


2u(k 2 -K 2 K 5 )Q 4 


K 2 (K 2 -kQ 4 ) 
2uQ 4 K 5 (k 2 -K 2 K 5 ) 


Q x f SCKj-kQ^ft-Lj) 

2K7 1 ^ 


Q, 


+ (Kj - Kg) < - uK + 2 k_ ( t_L l) 


K4U 


q; 


Q 2 (k 2 - K 4 Kg) 


■} 


L,K 4 f 3(t-Lj) u(K 4 -HsQ 4 n Pj 


2K 5 Q 4 


{- 


Q 4 u 


Q„ 


(k 2 -K 2 K 5 ) 


2 K n 


3u (t- L 1 ) + — ^ 


2 /t ,2 


(k + K 5 Q 4 ) k (K 4 - kQ 4 ) 


K„ 


K 5« 4< k2 - K 4 K 5> 


} 
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Kr, 


K. (K^ - kQ„) 

L 4 = (Kj - K,) i 4~- — (L, + Qo) + — 2 ° ~3T 

4 1 ^ u Q 4 (k 2 - K^K 5 ) ^4 3 3 Q 4 u(k 2 - K^K^) 


L 7 (2k 2 - K 5 kQ 4 - K 2 K 5 ) 


4 5' 


sin i cos (K^-f) K 4*^l 

^2 / o 2 • + 2 

sm i sirr (K^-f) uQ 4 


1A 


L„ = 


Q 1 K 4 u , ^ , (K 4 - kQ 4> K 4 

<4 oc 2 - 4 k 5 > ‘ 008 1 <l 3 + ^ ' Pl 


(K 2 + kQ,)u 


+ Lrr ~n 9- “ Qo 

7 Q 4 (k^ - k 4 K 5 ) 2 / “j” 

4 4 & VI - sin i 


sin i cos (K^-f) 


sm (K^-f) 


The ordinary differential equations will be developed here by taking the time differentials 
of the transf or mations and substituting in the original ordinary differential equations of r, 0, 
<|>, u, v, w, m, X, p, and X_. One notes that all terms will contain F as a multiple. These 
ordinary differential equations are presented for reference below and are used in the later 
discussions: 


w 


F h 

m~ A ( xj 

F r 2 X. 2 
m i(x) 


v 2 

r 3 ” 


w 2 sec 2 0 


w 

2~ 


sec 0 tan 0 


F r 2 Xg cos 2 0 
m A ( x ) 


(22) 


r 

0 


v 

“2" 


w 

cos 2 0 


= -X 


2v 

2w 


tJV X 

l^- p 2v- 2 


= _X 1 ^ 


sec" 5 0 + X£ 


2w sec 2 etan e 


p 3 


sec 2 e 


11 



s fs (v 2 ± sec 2 0) . 2|j. , 2 to? sec 2 , tan , 


2v 


+ p 2 "3 + p 3 T sec 6 " m 


2 w 2 Q F r(X^ + x2 C os 2 e) 
A(x) ' 


2 ^ 

P 2 = - Xj 2 sec 2 e tan e + x 2 ^ ( 2 sec2 6 tan2 6 + sec;4 6 ) 


K 3 

m 


9 „ r 2 sin 6 cos 6 

Wo 2 0 , Q F O 

" p 3 T 2 sec 6 tan 6 + — ApT) 

r 


A(X) 


m 


The ordinary differential equations for K.(K.) are easily derived by combining Eqs. (20) 
and (22) in the following order: Kg, Kg, K^, K ; - , Kg, K^, and K^. One then obtains the 

following expressions: 


Or - r 
*2 ~ T " 
m 


A(X) 


FQ4Q0 cos 2 e 

K 3 = m A lx)" 

* F + 

h ~ "m ^(5cpq 

F (^U+QgV+KgQg 
K 5 " ■ 2 “m ST(x) 

F tane q| /KgQg -vcos 2 e Q g 
K 6 = m A W \ (K 2 - ) 


F J Q 4 K 3 tanS (Q 2 - v cos 2 6 Qg 

*7 = mA(X) 


{ 


K 4 (K 2 -K 2 ) 

k| (K 2 - kQ 4 ) Qj - (QgV + Q3K3) Q 2 u (K^ + kQ 4 ) 

Q4 (k 2 -k|k 5 )k 4 


} 


(23) 
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where 


sin i cos (K^ - f) 
v = — — 

1 - sin 2 i sin^ (K^ -f) 

K 2 

2 _ K , 2k ^4 

" «r 


cos 6 = ^ 1 - sin 2 i sin 2 (K^ - f) 
sin e = - sin i sin (K^ -f) 

K, 


(24) 


% = ^jT W + ml~U) |^ P 1 Q 1 + ^Sa “ Q 3 L 6 + 

2K« 2 2 2 2 

+ Q 1 Q 3 -q^. - Q2Q3 2K 3 tan e " uQ 4 (Qg + Q 3 cos 0) + cos 0 sin 0 vQ 3 

where and A (x) will be defined later. All other variables are K's, Q’s, and t. 

To obtain the L’s, one observes that 

t - 9S(K, Q,t) 
l 9K. 


Therefore, the general equation for the time derivatives is given by 


L. 


3 2 S 


1 aiq 3K 


Q + Js 


QK^dQ 


aic at 


(25) 


where 


9K.9q 
1 


^m = X 


The 


j=i 

. 

9K.3K. ’ 9K.9Q. 

1 ] 1 ] 


9 2 S 

3K.9q. '‘j 
x 




2 

9 S 


and 


3 2 S 

3K.3t 

1 
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will 1)0 given in a document (previously cited*) which will be published in the near future. 
Remembering that the terms must be a multiple of F, we may write: 


P 1 


where 


Ql 


where 


Q = 


o 

0 

0 

0 

F r 2 sin 0 cos 0 
m 
0 
0 


w 


may be defined from the Hamiltonian, so that 


n as 

p l = “ff 


[ K 1 _K 2 +Q 1^‘ 
+ Q 5 ^ Q 6 K 3f^] 


K *4 

of 


k > 


Qg ^2“ sec2e tane 

q 4 


Q 6 _ " L 6 

Qrj = L2 cr -O- (t - Lj) 

, Qg, Q 3 , and Q 5 are obtained from Eqs. (21): 

c 

K 4 L 7 (k| - kQ 4 ) 

Q * =2uL *-ira' 


Q, 


4 5' 


2u 2 Q. 


k: 


+ ■ , i COS 


2k; 


5/2 


_! k-K 5 Q 4 k(kK| + K 5 I^Q 4 -2k 2 Q 4 ) 

k| (k 2 -K^K^) 


" 

= A' 1 

G 7 

% 


G 3 

_ Q 3 _ 


_ G 4 _ 


k2 " K 5 K l 


= A' 1 g 


A and g are given below. 


( 26 ) 


(27) 


(28) 


(29) 


( 30 ) 


*See footnote on page 3 . 
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2 

sin i cos i sin (K^ - f ) 

[ 1 - sin 2 i sin 11 (K ? -f)] a/2 

cos icos(K^ -f) 
sin i [ 1 - sin 2 isin^K^ -f)]' 

~ sinicosCKrj -f) 

1 -sin 2 isin 2 (K 7 -f) 


+ sin i cos (K^ -f) 

/l - sin 2 i sin 2 (K^ - f ) 

-cos i sin (K^ -f) 

- K \J 1 -sin 2 i sin 2 (K 7 -f) 


o 

cos isin(K 7 -f) 

-sin 2 isin 2 (K 7 -f) 

cos(K 7 -f)[ 1 - sin 4 i sin 2 (K 7 -f)] 
sin i [ 1 - sin 2 i sin 2 (K^ - f)f*^ 2 

- cos i sin(K 7 - f) cos(K 7 -f) 

K 4 [l - sin 2 i sin 2 (K 7 -f)] 


cos i sin icos(K 7 -f) 
j X - sin 2 i sin 2 (K 7 -f) 

3 

-K^sinicos isin(K 7 -f) 
[1 - sin 2 i sin 2 (K ? - f)] 3/2 

- sin 2 icos 2 (K 7 -f) 

[1 - sin 2 i sin 2 (K 7 -f)] 


- sinicos(K 7 -f) 

f o 9 

1 - sin i sin (K ? -f) 


2 

sin i cos isin(K 7 ~f) 
[1 - sin 2 i sin 2 (K 7 ^f) ] ^ 


-cosi 

[1 - sin 2 i sin 2 (K 7 -f)] 


L 7 + L 6 


l 3 + L g 


1 - sin 2 1 sin (K^ - f) 

cos(K 7 -f) sin(K^ -f) 
K^[l - sin 2 isin 2 (K^ -f)] 


Q x uK 2 


k 4 (K^ -kQ 4 ) 

1 - Ln cos i + Pi it n 

4 K 5 ) 3 1 Q 4 (k 2 -K^K 5 ) 


u(K^ + kQ-) 


I 


is obtained by iterating a solution of 


t - L, 


Q 4 u k 

-^- + ^ cos 


■1 K-KsQi 


" 2 - K 4 K 5 


with t given. 

a (x ) is defined as 

A ( x ) = \/\ + ^ r 2, + x| r 2 cos^ 0 (32) 

where Xj = Q 1 is given by Eq. (29) and and X 3 are and Q 3 given by Eqs. (30) and (31). 
The new Hamiltonian is given by 

” = -% ♦ <33 > 


The Delaunay variables may be given by 
o 2 

a . 2k _ 2 „ °6 
a 5 ~ r u Z2" 


? 2 2 2 
= v + arj sec 6 


arj = w 


i a™ tan 0 

ao = $ - sin” 

/~2 2 “ 


= cos 


V a 6 “ 

2 i 
a 6 ~ kr 

r J fc2 ‘ “5 ff 6 


o' sin 0 


f a 6 “ a 7 


, „ ruk -1 k -V 

* + B i " - T72 + 008 7f=z 

> 2 -Ve 


K we compare these equations with Eqs. (20) and the first equation of Eqs. (21), we see 
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( 35 ) 


K 5 = “5 
K 6 = “3 


K 7 = °4 
*3 = “7 


K 4 = a, 
*1 


6 3/2 
a 5 

*nr~ 


We let 3/2 

Or 


*2 “ a 2 -TT- ’ 

and then our generating function becomes 


a 3/2 q 3/2 

S = ~ a l -Tc- 4 + a 2 irir- Q 7 + < a l + a 2> 


u Q 4 « 5 

— Z 


1/2 


+ cos 


-1 k -*5 Q 4 


1 2 2 
k -o 5 a 6 


2k “6 


" a 7 q 3" \ / _a 5 + T2 ra - 

^4 ^4 


(36) 


+ Q 5 sin~* 


sin i sin ( a ^ - f) 


% 


sin i cos (ct^-f) 

^Jl - sin 2 i sin 2 ^" - f)~ 


where 


+ 


-a^ + sin 


cos i sin ( a 4 - f) 


I -sin 2 i sin 2 ( a 4 - f) 


/ 2 2 2 , A 
Qfn — k Q* 

sin i = v 6 .. 7 > cos f — -S ^ 




q 4v A 2 To^ 


, u 2 = 


- a 5 + 


2k 6 

Qa ’o!" 


We treat S as before, so that 


*3/2 uQ 1/2 k -« 5 Q 4 


(37) 


3/2 


^2 = 


r> 3 / 2 

Q 7 . a 5 


t + 


% = - % 
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*T-«! 


+ Q c 


cos i cos ( a £ - f) 
sin i jl - sin 2 i sin 2 (a^ - f)J 

sin(a 4 - f) cos(a 4 - f) 
ot g jl - sin 2 i sin 2 (a ^ - f)j 


3/2 “ Q 5 


- Q 3 


sin ( a 4 - f) cos i 


2 2 

oig /I - sin i sin (a^ - f) sin i 


( 38 ) 




a g sin i cos i sin ( a ^ - f) 


1 - 


+ Q, 


2 2 
sin i sin 


cos 1 


6 r . 2. . 2 

1 - sin l sm 


(“4 " f) ] 
(g 4 - f) ] 


372 + Q 5 


sin i cos - f) 


/ . 2. . 2 , ,, 
- sin i sm (< ot ^ - f) 


% = 


^4 a 6 (a 6 ~ kQ 4 } 

2 u Q 4 (k 2 - a 5 o!g) 

+ (a l + a 2 ) 


3 a 


1/2 


27 + Q 1 2 u + 2 


( Q 7 a. > 

-7- - r -t 

t 2 crk k 


uQ 4 q // 2 < 4 k + “ 5 « 6 Q 4 - 2 k 2 Q 4 > 

13^ + ' 2 «j /S q 4 u(?:« 5 4) 


3 a 


1/2 


‘5 Q 7 _ t , . „ « 6 a 5 u(a 6 + kQ 4 ) 

/3 p- s (®o - ®i i. ) ft A o o o 

5 2 2(7 lk 4 2 a 2 Q 4 (k 2 - a/ 6 ) 


+ Pi 


q 4 

g 6 (g 6 " kQ 4> 

+ Q x 

i 

2 

^ CD 
53 

H 

0 

i 

1 

g 5 

2 a 5 Q 4 (k 2 - a 5 a 2 6 ) 

2g 5 2 Q 2 a 5 (k 2 - a 5 a 2 ) 


^6 " 


(a 1 + a 2 ) a 3 5 /2 a 6 (a 2 6 - kQ 4 ) 

2^ 2 
u Q 4 k(k - « 5 « 6 ) 


a. 


+ Q 


1 uQ 2 


-<v 


sin i cos(a^ - f) 


^1 - sin^i sin^ (a^ - f) 


- (/3 7 + Q 3 ) cos i + 0 4 


2 k2 - g 5 a 6 ~ g 5 kQ 4 
uQ 4 (k 2 - Qfgffg) 
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0 R = Q 


a 6 U 


„ a 6 (a 6- kQ 4> 

1 „2 2 2 , „ ,.2 2 > P 1 " ^2 / = s 

Q 4 (k - a 5 a 6) Q 4 a 5 a 6^ ^/l- sin 2 i sin 2 (a 4 -f) 


sin i cos - f) 


u (a* +kQ 4 ) 


- 03 7 + Q ) cosi + ^ 4 — -g *- 

Q 4 (k -a 5 a 6 ) 


a l = 


' a 2 + 


k 

1372 


5 L 


Qi 


2 

ff 6 ” kQ 4 

Q 3 


+ P l u + 'T2 _ ^4 
y 4 


We repeat the operations for determining the ordinary differential equations. Eqs. (22) 
are still applicable to this problem and are used: 


a l 


FQ 2 Q 3 cos 2 e 
mi (x) 


FQ^(vQ 2 

a 6 ~~ mi (x ) ffg 


a 5 = " 2 


F (QjU + Q 2 v + a? 7 Q 3 ) 
mi (x) 


a 4 = 


2 2 
F Q 4 tan 0 (« 7 Q 2 “ V Q 3 cos e ) “7 

“2 27 


mA (x) (a 6 - a ? )a 6 

a 2 (a 2 -kQ 4 ) Q x - (a 2 + kQ 4 ) Q 2 u(vQ 2 + « 7 Q 3 )] F 
a Q Q 4 (k 2 -«6 “5) mA 


a 3 = 


F tan e Q 2 (a^ Q 2 - V Q 2 cos2 e ) 
mA(x) (a 2 -a 2 ) 


(39) 


F a i ( x ) k 3a 2 

a 2 = n 2 3/2 "2^7 a 5 

Q 7 « 5 5 
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Q 


T " A(X) + SMx) 


P 1 Q 1 + Q 5S " Q 3^3 


2 a r 


+ Q l Q 2 qJ +Q 1 Q 3 of -Q2 Q 3 2 “7 tane ' uQ 4 (Q 2 +Q 3 COS 0) 


where 


+ cos0 sin0vQ„ 


3 

^ “ 5 


v = 


sin i cos (a^ - f) 
^/l - sin^ i sin^ (a^ - f) 


2 2k ^6 

u ~~ + O ”2 

5 Q 4 Q* 


/ 2 2 

cose = /I - sin isin (a^ - f) 


sinO = -sin isin (< a ^ - f) 


The terms p^ and A(x) will be defined later. All other variables are Q f s, a' s, and t. We 
use the same procedures as before to determine /3 r s; that is: 


£4 


as (a, q, t) 


9 a. 
1 


so that 


A = 


9 s a a s 


da. da 

X 


a + 


9a.9Q 


Q + 


9 2 S 

8a. at 
1 


(40) 


00 o 

S g 9^S 

Again -r — jr — , - — , and -r — kt will be given in the report cited in the footnote on page 3. 
da. da. 9 da. 9Q. ’ 9of.9t ^ 

i] 13 J_ 

Equation (26) is used to define Q and p^ is given by 

"N 


3/2 


f ots 


p l = < l-^l +a 2 ) — + Q ll^3 ^- sec 0tan<t> + Q 5 J +/3 3 a 7 ^ 


2 

sec 0 1 


q: 
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Q-7 " 


3/2 


Q x = 2 u /3g - 


3 uai 


1/2 


f q 7 

a 9 a^t 

l ^ O' 1 


(«1 + 0f 2 ) a 5 /2 

kQ, 


2 r 


2 a e(‘V kQ 4 ) 

T 2 2 — 

k -a 5 a 6 


+ a 6 /3 4 (a 6~ kQ 4 ) 

Q 4 (k 2 - «g «g) 


We use the same formulations as before for computing and A (x) where 

sin i cos (a 4 -f) 

“ n. 


QTg sin i cos^ i sin (a^ - f) 


1 - sin^ i sin^ (< ot 4 - f) 


cos i cos (a^-f) 


W 2 ~ 


i jjl — sin^ i sin^ (a^ - f)j 

- sin i cos ( a 4 - f) 

2 2 

1 - sin i sin (a^ - f) 




- sin 2 i sin 2 (a^ - f) 

- cos i sin (or^ - f) 

/ 2 2 

ffg sin iy 1 - sin i sin (a^ - f) 


-1 


(41) 


(42) 


2 

cos i sin (a^ - f) 

2. . . , 
cos isinicos(a 4 - 

f) 

- sin i cos (a^ - f) 

y/ a l ~ [/■ ~ sin2 i sin 2 (a^ - f) 

/ 2 2 
1 - sin i sin (a^ - f) 

/ 2 2 
+y 1 - sin i sin (a^ - f) 

cos ( a ^ - f) [l - sin^ i sin 2 ( at ^ - f) j 

- ffg sin i cos'* i sin (a 

4- f ) 

2 

ag sin i cos i sin (a^ - f) 

T 2 2 I 3 / 2 

sin i 1- sin i sin (a^ - f) 

jl - sin^ i sin^(of^ - f)J 

ttt 

r 2 2 i 3 / 2 

1 - sin i sin (a 4 - f) 

- cos i sin (a^ - f) cos (a^ - f) 

2 2 

- sin i cos (a^ - f) 


- cosi 

2 2 

ffg 1 - sin i sin (a^ - f) 

£l - sin^ i sin^ (a^ - f) 

1 - sin^ i sin^ (a^ - f) 
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g = 


Pa + P‘1 


Pn + Pr 


COS 1 


j^l - sin 2 i sin 2 (a^ - f)j 

sin (a^ - f) cos (a^ - f) 
|l - sin 2 i sin 2 (a^ - f) 


Pa ~ 


Q x ag u 

Q 4 ( k2 - “5 4 


+ firj cos i + 


a 6 ^6 ~ ^Q 4 ) 
Q 4 (k 2 - « 5 a 2 ) 


u (cifg + kQ 4 ) 

Q 4 (k - a ^ flg) 


2 2 

For the Poincar<§ variables, we wish to eliminate the singularities of K A = K Q 
2 2 *± o 

(zero inclination) and k n K^K 4 (zero eccentricity). If we view Eqs. (13) through (17) and 

Eq. (19), we see that this cannot be done with the functions used. A look at Eqs. (15) and (16) 

suggests adding or subtracting Kg and to try to eliminate the K 2 = Kg singularity. 

When this is done, one has 


Kg - = 4> - cos 


-1 K 4 ' 


r /k 2 - K 2 K 5 


+ sin 


-1 v sin 9 


K 4 + K 3 


(43) 


where >0 and Kg > 0 by convention. 


3/2 


5 k 

We next note that if we choose K^ = - a ^ ^ — , we remove the factor from 

Ctr- 


cos 


-1 k - K 5«4 


k 2 - k 5 k 2 


term and have possibilities of combining it with the K 4 - Kg term. Also, 


note that the /3j term that results is a position- velocity term (no Lagrange multiplier -- part 

of the unadjoined two-body problem). Assume again an addition or subtraction of the /3 1 term 

2 2 ^ 

with the K^ - Kg term. This results in the elimination of the k = K^K 4 singularity. When 
this is done, one has 


as 

da* 


= <{> + sin" 


■1 vsin 0 


k; 


3/2 


jrl/2 

ru K 5 


k 4 + k 3 


+ cos 


Ti — is! — 1. 

L b- J. tA-I j 


k+K- -K 4 


(44) 
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We have now established one new variable. If we solve Eq. (15) for tan (- tan e), 

we may write: ^5 


sin K _ v sin 4> - w tan e cos <\> 
6 


/ 


k? - k: 


cos K = v CQS 4> + w tan 9 sin <|> 
6 


/ 


k \ - k: 


When multiplying both sides byy'K^ - w, we obtain: 




/K, - w sin K, = v sin <t> - w tan e cosjj = 

V * o /-fp ; 6 

y K 4 - w cos Kg = 


k 4 + W 
v cos cf) + w tan 0 sin 


(45) 


= a. 


y K 4 + w 

If we operate on Eq. (16) in a like manner, we obtain: 


«6 =/* - K 4 K 5 /2 sin (K 6 - *9 

(K 4 cos 0 + K 3 sec e) [ sin <f> (K^ - kr) - cos <j> K 4 ru] + v sin e (cos 4> (K 4 - kr) + sin 4> K 4 ur) 


r (K 4 + Kgj/k + K 4 K ^ /2 


(46) 


=/ 


k-K 4 K^ /2 cos(K 6 -K 7 ) 


[(K 4 COS e + Kg sec e) (cos 4> (k| - kr) + sin <f> K 4 ru) - v sin 0 (sin <j> (K 4 - kr) - cos 4> K 4 ur) 

r(K 4 + K 3 )/k + K 4 K | /2 

Let n 5 = Kg, and we have given a 5 , a 3 , a 4 , a g , and from Kg, K 4 , K 3 , Kg, K ? or rather 
(u, v, w, r 0, 4>) where the a's have no zero in the denominator at zero inclination or 
eccentricity. 
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Our transformations become: 

2 , 2k sec^ 0 

«5 = ~ u + — - j — 


^3 = 


v sin <\> - w tan 0 cos <t> 


2 2 2 
v + w sec 0 + w 


a 4 = 


v cos 4> + w tan 0 sin 4> 



2 2 2 
v + w sec 0 +w 


K a =V v + w sec 0 


a 6 = 


2 

(K^ cos 0 + w sec 0) (sin $ (K^ - kr) - cos <j> K^ur) 


r ®4 * %/k + K 4 o| 


/ 2 


+ v sine [(cos 4> (k| - kr) + sin 4> K^ur) j 
r(K 4 + w)^/k + K 4 a^ 2 


“7 = 


(K^ cos 0 + w sec 0) (cos 4> (K^ - kr) + sin 4> K^ru) 


(K 4 + w)^k + K 4 4 /2 


v sin e[sin <j> (k| - kr) - cos <t> K 4 urj 


r (K 4 + w) J k + K 4 


1/2 


Next we write the I^'s as functions of a i (i = 3, 4, 5, 6, 7) 



K, 


.22 

k - a 6~ a 7 

1/2 


2 2 
a 3 " a 4 


K. = 


i 2 2 

k ~ a 6 - g 7 

/v 1/2 


K 5 = a 5 


( 48 ) 


-1 a 3 
K c = tan 1 -4 


K^, = tan' 


a, a™ - a. a 


-1 “3 7 


4 6 


a 4 a 7 + a 3 a 6 


We next note that, if S (K, (5, t) is a solution of the Hamilton- Jacobi equation, then so is 
S (K, 73, t) + S* (K) .Let 

.3/2 


K 1 g l “TT 


, 3/2 


K 2 “ a 2 "TT 


in S and note S = S (old) + S 


as 

ao?i 


3/2 


t - 


“V / 2 -1 

4- + cos 


k - a 5 Q 4 


/<4 + “?) (2k _ “ a 7 ) 


8S* 

‘cTaT 


||l-= K 6 -K 7 = tan - * — - tan -1 “ 3 

67 “4 « 4 « 7 + «3 a 6 


= tan 


-1 


Clc 


so that 


S* = (a 1 + a 9 ) tan 


•1 “6 
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The resulting new S is now written as S^ new ^ = + S*. 


v 3/2 


S =■ 


Q, 


(a 2— ~ “1*) + ( a l + a 2 ] \ 


1/2 

UQ 4°5 + tan - * 

E arj 


+ cos 


-1 


k - 0:5 Q 4 \ / k - ac - a? 


y^ + a^) (2k - a* - a% 


-Qa 


6 " 7 


4 /2 


a 3 - “4 


- 4 T 


'" 5+ ^ 


-1 


a 3 


- Qe tan cn 


( 49 ) 


K 4 sin i cos (K_ - f) _i 

- Qp . i + Qc sin [ sin i sin (K„ - f)] 

Vl - sin 2 i sin 2 (K ? - f) 

cos i sin (K^ - f) 


+ Qg sin 


vl - sin^ i sin^ (K 7 - f) 


where 


sm l = 



cos f = 


[ k ~ g 6 - °^\ " a 5 k ^ 


K 4 = 


a 5 Q 4 y^e + a 7^ ^ 2k _ a 6 " 

. 2 2 
k - a 6 - °7 


a 


1/2 


( 50 ) 
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, „ g 3 g 7 - g 4 a 6 

7 " a 4 a 7 + a 3 °q 






vsin 0 

i 2 2 

k - a 6 - a 7 2 2 

2 1/2 a 3 - “4 

g 5 


„ 1/2 

Q 4 u « 5 ' 


-1 , 

+ cos 1 2 2 

2k- aig - 


3/2 1 

“5 Q 7 . „ 

__ — + t +/3 1 



a 4 +a 3 


/ / , 2 2 

/2 2 N ( 0 k_a 6 g 7 2 2 

a 3 + ® 4 ^ 1 2 1/2 a 3 ~ a 4 

\ a 5 


cx 3+ a 4 




- Q x (k- a 2 - a 2 ) u Pj k (k- a 2 - a 2 ) [(k-a 2 - a 2 ) - ka 5 qJ 


„2 2 , 2 2 , 2 2 , 
Q 4 a 5 (“ 6 + « 7 ) (2k- a 6 - atrj) 


2 „ 2 2 W 2 2. 

« 5 Q 4 k(2k- a 6 - a ? ) (ccg + a 7 ) 


^3 9Z 1 , 9Z 1 

~ T /2 - 9K . 1/2 + 9i . 2 2 

a 5 4 a 5 k ~ a 6~ a 7 


a 3 +a 4 

1 2 2 

k -« 6 -« 7 


- g 3" a 4 


27 



■■mil 


dZ 

8(K„ -f) "372 


u [(k - <x\ - a 2 ) 2 + kQ 4 « 5 J 


7 " I! a*' * Q 4 (a 2 + a 2 ) (2k - a 2 - a 2 ) 



a 6 

/ 

+ < 

. 

“ 7 . 



ku 


kQ 


! 1 [fr ~ g 6 ' " 


5/2 , 2 2 X „3 
“ 5 (a 6 + a 7 )Q 4 


) 2 - 5 kQ 4 ] 9Z 

+ 9(K ? -f) 


k (k - a\ - a 2 ) - a\ Q 2 


2 ~2 , 2 2 , 

“5 Q 4 (a 6 + a 7 ) 


=2 


(“l + a 2 ) a \ /2 (k ~ a l ~A ] 


(k _ a 2) 2 


-kQ, 


uQ 4 k (a 2 + a 2 ) (2k - a 2 - a 2 ) 


\-OLc 


Q 3 

Qj (k - a 2 - a 2 ) 


8 Z 

1 2 2 


1 8Z 

8i 

1 a 3 + a 4 


1/2 

a 5 

„2 

u Q 4 « 5 

~ T 72 9K, + , 
a ^ 4 k 

2 2 1 

- a 6- a 7 

-If "V 

, 2 2 

0 k “ a 6 ” a 7 2 

—172 a 3 

a 5 

2 

- a 4 


9Z k(k “ a 5 Q 4 } + (a 6 + a 7 } (2k ” 4 ~ A ] 

2 2 


^ K 7 ‘ f) ~ u l z{ 2 u Q 4 («| + a 2 ) (2k - atg - a 7 ) 


a . 


2 2 \ a l +a 2 " 9(K„ - f) 

ffg + a ? \ v 7 '> 


9Z 


, 1/2 


n 7 


- a £ 


SCtfc / Qr 7 \ u 

% = -air - ( a 2 — " "i 4 ) " Q i + 


% a 3 + ^4 a 4 


Q 4 + 

P 1 *i + 


Q, r k - at - a* o ? 1 sin 1 cos - i) vs. 


sin i cos (K^ - f) (k - a?g - a 

sin 2 (K 7 - f) 
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I 


a 3a 5 ( S .\ 

VlF “27-“!* + 


«5 ^ 2 Q 4 («! + « 2 ) Q 3 (k - a 2 - a 2 ) «g Q 4 - (k - «g - a 3 ) 


8 Z 1 
8 i 2 a c 


2 2 
a 3+ a 4 


2«s /2 2 u « 2 5 Q4 

9 Z k ' a 6" a 7 9 Z k ' a 6‘ a 7 


k _ ff 2 „ 2 8 K. „ 7372 “ 

, k 6 7 2 2 4 2a 5 

2 172 a 3 " a 4 

a 5 


C„ - f) 3/2 „ 

7 ' a 5 ' uQ 4 


where 


92 _ Q 

8 i _Q 2 


( 2 2 \ 2 2 2 

k - ag - jcos (<)> - Kg) cos i + sin i sin ( 4 > - Kg) 


cos i 


- Qg sin (<j> - Kg) + Qg tan i cos (<|) - Kg) sin (<(> - Kg) 


„ 2 2, 

9 Z (k-a 6 -a 7 ) 

b(K 7 -f)' = — 172 tan 1 sin 

5 


W-k 6 ) [ 


2 2 2 
cos i + sin i sin (<(5 - K. 


+ Q 5 Sin i cos ( 4 > - K^) + Q g 


2 2 2 
cos i + sin i sin ( 4 > - K^) 


9 ^ - = -Qg sin i cos ( 4 > - K 0 ). 

One notes that, when a Q = a A =0 and a R = a „ = 0 , there are singularities in the above equa- 
tions. 6 4 57 


If one defines the following: 


“3 £3 + a 4 P 4 2 2, ^ 9 Z 

A 1 = 2 = (a 3 + “4 } Q 3 + 8 r 


/ 2 2 
/ “3 + 

' , 2 2 

k ” a 6 ” a 7 2 2 

2 1/2 - a 3 ~ a 4 

a 5 


^2 a 4 @3 “3 ^4 "^6 + 9 (K 7 - f) 


A 3 a 7 / 3 6 ' a 6 /3 7 a l + “2 ‘ 9 (K ? - f) 

Pl ku Qjk [(k - a 2 - a 2 ) 2 - c*gkQ 4 ] 

A 3 = ^372 + 4 

5 a 5 ^4 


, /. 2 2, 2 „2 
9 Z k(k-a 6 -a 7 )-a 5 Q 4 

„ 2 Q 2 


l 




, , 1/2 , 2 2 X 
(o' 1 + Q' 2 )a 5 / (k-a 6 -a 7 ) 

2 2. 2 
(k-Ug- » 7 ) 

L “5 

uQ 4 k(2k- ag 

i ' 
;-« 7 ) 




2 2 , 

(k-a 6 -Q! 7 ) u 

2 „2 /01 2 27 

a 5 Q 4 (2k-a' 6 -a 7 ) 


(k-a^-Q^) [^(k-a 6 -a 7 ) -ka 5 Q, 
agQ 4 (2k-ag-a 7 ) 



and uses these along with the equations for p-, p„, and p g , then the system is non-singular 
at a 6 = a i = ° and “3 = a 4 = °* 
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a 1 / 2 [v + QgtKg + ^cos 2 ^] 




2 2 
a 3 ' a 4 


mA( x ) K 4 (K 4 + w) 


/— 

cos <f> 

sin 4) " 


/— 

(vQ„ + wQo) Q^u + Q t 

< 

.sin <J> 

cos <J>_ 


(K A cos 0 + w sec 0) 

4 0 4 14 

v_ 


+ w 



2Q 4 K 4 (vQ 2 + wQg) 



K 4 uQ 4 

( 

1 

Sf 

+ 

oq ^ 

W 

1 

* * 

+ Q 4 cos e (vQ 2 + (w + K 4 )Q 3 ) 


+ sin e K 4 Q 4 Q 2 



K 4 - kQ 4 

V 

1 

9f 

W 

. • t 


+ V 


-2Q^K^ (VQ 2 + wQg) 


(vQ 2 + wQ 3 ) Q^u + Qj 



VQ 2 + Qg (w + cos 0 ) 


[ (q 2 v + Q 3 w) (Q^a 5 - K 4 ) - Q x uK^ (K 4 + w) 

+ 0 1/2 7, “ T 727 

2 (k + K 4 a g ' ) 


= -2 


F(QjU + Q 2 v + Q 3 w) 
mA( x ) 

3ca„ 


_ kF o- A(x) 

3/2 2 " 2 a- a 5 

1X1 y 


( 54 ) 
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3 (o' 1 +« 2 ) _ Fk 

1 = " a 2— 2^— “5-.372^7 x) 


P 1 Q 1 + SSz + Q 3 Q 6 


+ ^~ + Q l Q 3 ^ L - 2 S Q 3 wtan0 - uQ 4 (Q l + Q 3 cos2e) 


+ COS0 sinGyi 




FQ^cos (K^ -f) 


J t mA (x)K 4 


COS1 


iQ 2 


/ . 2 . . 2 ,„ 

' 1 - sin i sin (K^ -f) 


- Qg sin i cos (K^ -f) 


cosisin (K 7 -f) sin (4>-K 6 ) 

[ K 7 sin i cos (K^-f) \ ~~ sin i cos^-K^) *t 


(55) 


The following equations may be used to evaluate the Q's needed in the problem and give 
information for A (x). 

Adding A ^ and A^ of Eqs. (53), one obtains: 

Q 6 = W ^ ‘ A S = a l +a 2 +a 3 P 4 + V? " «7 P 6 ‘ «4 P 3 
From equations for and of Eqs. (53), one obtains: 


g^ = sin i cos (4 -Kg) + 




k : \ .£.% 
„ 1/2 3 4 


Sta(*-K 6 )[ V Q 6 &^] 


9 9 [cos^i+sin 2 isin 2 (4>-K fi )l ( a o+^) 

= (« 3 +^) sinicos(<j>-K 6 )Q 3 - 


V ,3 
8Z 




9Z 9Z 

Next we eliminate — yy and -gpfrom A^ and by Eqs. (53) for A^ and and 

obtain: ^ 
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, 2 2,2 

i/ 2 , 2 2 , ( k - a 6 ~ a 7 ) 


g 2 A 4 " 


(a l +a 2 )a 5 (k-« 6 -« 7 )^ ^ ~-kQ 4 

uQ 4 k (2k - a 2 - a 2 ) 


* k (k - a 5 Q 4 ) + («g + a 2 ) (2k - a\ - a 2 ) («g + <* 2 ) 


«5 ^ 2 uQ 4 (2k - a 2 - a 2 ) 




, 2 2 , 
(a g + a 7 ) 


k - a e " a 7 sin 1 cos ^ _ V cos i 

n 2 + 173 1 

uQ 4 a 5 « 5 ' a 5 ' 


„ 0 “S 72 , «7 2 “5 /2 «4'“l + “2» 

g 3 " 2a 5 P 5 ” k ^ a 2 a- a l^ " ku 


2(k-a 2 -a 2 ) * a 2 Q 2 - (k - a 2 - a 2 ) 2 

+ 1 / 2 * *»-*!-— q 2 0 . 

a 5 uQ 4 ua 5 Q 4 


„ (k- a 2 - a 2 ) sinicos (0 - Kg) ^(k-a 2 -a 2 2 g 

+ ^ ' 71/2 +Q 3 ~T72 a 3' a 4 

a 5 \ a 5 


A 3 = a l + a 2 * A 3 * 


Next we write: 


’il M k 


is a 3 x 3 matrix. 
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or 



g i 


V 

[B] _1 

g 2 

= 



g 3 


_ Q 3 


and a 2 (x) = Q 2 + Q 2 Q^ +Q4Cos 2 9q| . 


-[cos 2 i + sin 2 i sin 2 (4) - K g )j (a 2 + a 2 ) 


- (/ + ah (k - a 2 - a 2 ) (a 2 +a 2 ) 

7 6 7 - 6 f7a 7 - sin i cos (<j> - Kg) 

5 


uQ 4 flf 5 


2 _2 2 2, 2 

a 5 Q 4 -(k-a 6 -a 7 ) 


uQ 4 a 5 


2 2 

(k-ag-a^) sin i cos (<{> - Kg) 

1/2 


(a; 2 + o^) sinicos (<}>>Kg) 


2 2 
a 6 + *7 

-T72- ( 


k -“6- a 7 


T7T 


a 2 a 2 




-COS 1 
“12 2 
*3 +a 4 


“2 2 “ 
<V a 4 


-uK 


A 


K +a 7 ] % 


sin i cos (<f>- K g ) (a 5 Q^-K^) 

X 


~^ /T K + a i ) 


u 


sin i cos (<J>~Kg) 


Q 2 a 

W 4 5 


sin i cos (<j>- K g ) [cos 2 i + sin 2 i sin 2 (4> - Kg)] (a^Q^-K^) K 4 [cos 2 i + sin 2 i sin 2 (<f> - K g 


1/2 Jl, 2 2 , 

cos i a 5 7 Q 4 (Qfg+ap 


cosiQ 2 ^ 5 
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The equations for p n , and Q are given here for reference: 
-L 5 


“5 3a l /2 Q 7 

P 1 " Q7 ^5 + 2kQ. ^l 1 ■ a 2 ~ ) + Q 1 "2QT 


^3 a 3 + ^4 a 4_ %_ 

2Q 4 


4Q„ 


k ■ “'e ■ a 7 2 2 

„T72 ' “3 ' “4 


. . 1/2 /k ~ a 6 ~ a 7 2 2 

cin i /v 7 I . _ _ /v _ rv 


2 _ 2 . 1/2 


Q 2 sin i dfg' “ I 175— L ' a 3 - °ij (“4 cos t> + <z 3 sin 4 > ) (k - a & - aj a g 


2 /(Og+Q'^) 


2 2 , 


2 2 


L a 6 ” a 7 2 2, ... . . .. .... ^2 

Q^^/l - sin 2 i sin^ (K ? - f) 

1/2 ” “3 ” a ^) cos + «3 sin 9 ) 

°5 


q* = 


, 2 2 

k " a 6 " a 7 

Qo — ~T 7 o sin i sin ( 4> - K r ) + 


e 5 " ^3 ”172 

a 5 


2 2 

(k - a 6 - (Xrj ) (/3g cos <t> - )3 4 sin <j>) 

. 1/2 / k -“|-“7 2 2 

“5 r—jn V *4 



jAj sin ( <t> - Kg) + Qg cos (<t> - K g )j 


We may now proceed to determine A (x) and H as in the two former cases shown. Here 
we solve for <t> instead of (r) as before. 
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m. MODIFIED DELAUNAY AND POINCARE VARIABLES 


This section presents a method of introducing the three known constants of motion for 
the full problem as coordinates of the base problem. We will start with the Hamiltonian of 
section II: 


H 


/ 2 2 

x I v , w 
M 3 + 3 2 

\ r r cos 0 



v w 2 F 

+ p 2 + p 3 ~2 sec 6 ” X 7 <r+ m 


where 

A(x) = ^ + r 2 ^2 + r 2 cos 2 e Xg 

The first step is to transform the Hamiltonian and consider only the base problem. The 
transformation is: 


S = P 1 X 1 + P 7 m + P 2 X 2 + P 3 W + P 4 r + V + P 6 4>. 


(57) 


The results of this transformation give the following relationships: 


\ ■ Q 1 

-u 

= Pj = 8S/9Q 1 


\ = Q 2 

-V 

= P 2 = 9S/9Q 2 


r = q 4 

P 1 

= P 4 = 9S/9Q 4 


m = Q 7 


= P ? = 9S/8Q 7 

(58) 

3 

II 

cO 

CO 

x 3 

= P 3 = 9S/9Q 3 


0 =q 5 

p 2 

= Pg = 9S/9Q 5 


0 = Q 6 

p 3 

= Pg = 9S/9Qg 
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The Hamilton- Jacobi equation becomes: 


= 3S _ 


at 


as Q 

0 Q„ Q 1 




Q, 


Q 


3S as „ Q 3 se ° 2 Q 5 tan Q 5 as Q 3 sec2 Q 5 

- ^ + a§- 


as 


as 

a^ 


(59) 


9Q 4 0Qj 


q: 




8Q 5 Q 4 


We note that Q^, t, and do not appear, so that by separation of variables we have 


as 


as 


as 


at V * dQ aQ 


= K 


k 3 * 


(60) 


We have proceeded to this point as in section II. Here, we will introduce the known con- 
stants. These constants come from the following relationships in cartesian coordinates: 

X. x x + x x \ = 0 

This integrates to form: 

X. x x + x x X = M (61) 

where M is a constant vector of three components. The resulting equations in our coordi- 
nates are: 


= sin 0 A + cos 0 B 




cos 0 A - sin 0 B 


(62) 


where 


M 3 p 3 


A = X 3 v - w sec 0 + pg tan 0 
B = P 2 “ ^3 w tan 0 * 


It can be shown by Poisson brackets that using M^, Mg, and M« as new momenta is im- 
possible. However, we may use the following: 


„2 2 2 a 2 _2 

= Mj + Mg = A + B 


(63) 


M 3 ~ p 3 

One notes that v, w, Pg , <t>, p^). We, therefore, proceed with the fol- 

lowing differential equations of the characteristic strip: 
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( 64 ) 


9 S Q3 sec 2 Q 5 tan Q 5 


s/aQg 


r? 2 ^ 
+ Q3 sec Q 5 


as as 

9Q 4 8Qj 

We integrate in a manner similar to that used in section n to obtain: 


- Q 2 sec 2 Q g 


0S \ 2 = K + — — 

^l) 5 Q 4 'q 2 4 


We may now rewrite the Hamilton-Jacobi equation with the above constants in the form: 


2 (< 

( K 1 - K^r + Q x l - — - k 


9S_ as 2 

9Q 4 aQj r ^6^4 


,2 2 


as as 


Sj Q 3 sec Q 5 tan Q 5 - K3Q3 sec Q 5 + ^ 9^ 


From this separation of variables, we may write: 


K 6 K 4 = (Kj - ^r 2 + Q x [-T--*) + ^-ur 2 


K 4 K 6 = ^4 Sec2 Q 5 tan Q 5 - K 3 Q 3 S6c2 Q 5 + 
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From Eq. (63) we have: 



+ 


(" Ji~ “ ®2 Q 3 se ° 2 Q 5 + K 3 tan %) 
( ‘ ^ Q 3 tan Q 5) ' 


If one expands this expression and solves for 


as 

aQ 3 ’ 


one obtains: 


(73) 


as 

8Q 3 


-^2 


as 

9Q 2 


Q 3 se(j2 Q 5 


K 3 aQg tan ^5 


Q 3 tan Q 5 


(K^ - Q 3 ) 


'sM - Q 2) - (K 4 K 6 + 
^ " Q 3> 

Formally we may write: 


(74) 


,(%) 


r Q 3 _Q 2 Q 3 ^” sec2 Q 5 + K 3 ^ tan Q 5 + H' Q 3 tan Q 5 

= o/ 2 _2. 


(A -S> 


dQ„ 


f® 3 '•/ M 7 'A - «3> - (K 4 K 6 + K 3 Q 3> : 
J ~ 


(75) 


dQ Q 


One notes that the last expression on the right does not affect the Hamilton -Jacobi 
equation (59), so that we may augment our transfer function S by this quantity. 

One notes that + K? appears in the S function and does not affect the Hamiltonian. 
We let 6 

K 3 = m 7 = K 7 

and 

A* A - A 2 

and drop the aster isk(*). Our S function then becomes: 
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K, 

Zjt + — 


'ir^ cos 


-1 k ~ K 5 Q 4 

y^-^Ks 


- Qj u - Q 2 v + Kg Qg - sin 


-1 Q 3 tanQ 5 


4-4 


+ Kg cos 


+ K^ sin" 


•1 K|- kQ 4 

q 4 ^ - 4 % 


K 7 I V K 4 K 3 K 6„_ 

K 4v^- K »v^5 


-1 K 4 sin Q 5 

^4 - <1 


+ cos” 


+ Kg cos 


-1 K 3 K 4 " K 6 % 

' / 4 - 4 v 4 T 4 


Our K T s are now defined as 
K l=- H o 


K 2 =aX 7 


K 3 = p 3 

k| =v 2 + Q 2 sec 2 Q 5 

K 2 

K _2k k 4 2 

5 ' c 5 " 4 ‘" 




K 6 ~ - ' 2 4 sec2 Q 5 tan Q 5 ' p 3 Q 3 sec2 Q 5 " p 2 v) 

K 7 = (X 3 V " h Q 3 sec2 Q 5 + Kg tan Qg) 2 + (p 2 - Xg QgtanQg) 2 + K 2 
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We may now record the time derivatives of the K’s: 


± 2 =^ A(X) 

m 

Kg=0 

F(vQ 2 + Q 3 X 3 )r 2 
K 4 = mi(X) K 4 


K c = -2 


HaW ( v ^ 2 + S Q 3 + Q l u) 


(78) 


• -K 6 ‘ Fr 2 

K 6 ' -KJ- K 4 + 

K 7 = 0 


x^ Xg 2 Q 3 tan 6 - v sin e cos e X 2 - Kg Xg - p 2 j 


^1 ~ mA(X) ^" x l(|( vQ 2 + Q 3 V + p l) _X 2 (" ur *2 _2Q 3 X 3 tanQ 5 + ?2 ) 
+ Xg ( -XgCosQg (v sin Q 5 - ur cos Q 5 ) - p 3 )j + -^-A(x) . 


We next record the L T s: 

uQ 4 v 
L 1 " * + TC 7 


5 iP' 


Q 7 

Lj-lT’+t-Lj 


COS 


-1 k -K 5 Q 4 


J* 


fc 2 -K 4 K 5 


,3 = Q 6 - Sin 


1 Qo tan Qc i 

1 6 0 + cos 1 


k 4 k 6 + k 3 q 3 




(K 1 -K 2 )K 4 (K 2 -kQ 4 ) Q x K 4 K 4 KgK 4 QgtanQ 5 


" 4 ” u Q 4 (k^ - K 4 K g ) 


uQJ 


- 02 ^- 


V (K 4 - Qg) 


K c 


+ (k| Kg + Kg k Q 4 -2 k 2 ) Q 2 tan Qg 


uQ, (k 2 -K 2 K 5 ) 


V(K 4 - Qg) 


V A (K? - A - *3) - 2K 4 K 3 K 6 Q 3 - Q 3 


K 4 (K* - Q 2 ) 


( 79 ) 
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K 4 (K 2 -kQ 4 ) K 6 ( K 2 + kQ 4 )u Qg tan Q 5 


K 4 u ... . n _______ 

^ = Ql Q 2 (k* - K 2 K 5 )" P 1 Q 4 (k 2 - K 2 Kg) Q 4 (k 2 - K 2 K 5 ) ^ K 4 (K 2 - Q 2 ) 


v (K l + Qg tan 2 Qg) Kg Qg tan Qg v J(K * - Kp (K^ - Kg) Qg cos L ? 


Q2 K 4 (K 2 - Q 2 ) + K 4 (K 2 - Qg) 


r 2 t ^2x (r7 2 


K ? (K 2 - Qg) 


L c = 


Qi 


-3(t - Lj)(K 4 - kQ 4 ) 


K 4 u 


Q. 


Q 2 ( k 2 -K 4 K 5 ) 


-p ■ 


K 2 (K 4 -kQ 4 ) 

3 u(t - Lj ) - 2 Q 4 + -i-| 


Q 4 (k -K^ K g )_ 


+ K. K„ 
4 6 


(K 2 + kQ 4 ) u 3 (t - L x ) 


Q 4 (k 2 - K 2 K 5 ) " Q 2 


1 

2 K„ 


L c = 


(Ki - K 2 ) 

K c 


3 (t L) ® 4 t K 4«4 

2 (t - V _ i r + m 2 


kQ 4 ) 


2uQ 4 (k -K^K 5 ) 


> + 


Ql K 4 K 6 (K -kQ 4 ) 
25 2 uQ 4 (k 2 - K 2 Kg) 


L. = cos 

D 


-1 K 4- k °4 


Q, 


A - K 4 K 5 


-i K 4 sin Q 5 
sin — + cos 


.1 K 3 K 4 + K 6«3 


/ K 4 - «3 


■M, - «s -M - 4 


L„ = 


.! k; 0 3 * K 4 K 3 K 6 
sin ■ 


K 4 A 2 - K 3/ K ?- K 6 
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Note that these relationships were used: 


cos 


-1 K 3 K 4 + K 6 ^3 


COS 


JA -‘Ssi/Kt-K? 

-1 K 4 K 6 * K 3 Q 3 


= sin 


-1 


K 


JA - A JA -A 

For this problem 


= sin"'*' 


Q = 


o 

0 

Fr^g cos^0 


mA (x) 


0 

0 

0 

0 


Again, we may write 


3 2 S__ ~ + 9 2 S 


aKj 3K 


k + g2 s 
Q + 


where the 




COS Ln 


K, 


A - A 


K 4 yK 2 -K 2 cosL 7 


K. 


7 k 4 -q 1 


9 2 S 

3K.9q 

is obtained by FORMAC and K and O are given before. 

Next we will develop the inversions needed to define A (x): 


Ql — 2L^u - 


(Ki - 1^) 


Kc 


K 2 (K 2 - kQ 4 ) 

3u (t - L-.) - 2Q 4 + „ n 

1 4 Q 4 (k 2 ^K 2 K 5 ) 


(80) 


> + 


K 4 K 6 ( K 4 - k< 5 4 ) 

Q 4 (k 2 - K 2 Kg) 
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o = _.JL L , K 3 Q 3 ta " Q 5 . 2vL 5 , K 
2 K 4 4 (K 2 - Q 2 ) Q 2 


6 


uv (K 2 + kQ 4 ) q| tan Q g 

, 2 „2„ , + ,„2 ^ 


K 4 Q 4 (k <i -K 4 K 5 ) K 4 (K 4 - Qg) 


- V 


- K 3 2 - K 6> - 2 K 3 K 4 K 6 S ' i S + V(K 1- K 2> 


^ < K 4 * A ] 


^5 


2(t-L x ) (K^ + kQ 4 ) u 

^4 + (k 2 - K 2 K 5 ) 


V Q 2 


Q 3 sec 2 Q 5 tan Q 5 K 4 (K 4 Kg + K g Qg) v/ K 7 (K 4 * Q 3* ' 2 K 3 K 4 K 6 Q 3 * K 4 * K 6 + 

- n n “ — + ” ' 5 9 


v (K 2 - Q 2 ) 


< K 4 - S'* 


From our equations for L f s, we have: 


sin = 


K 7 Q 3 + K 4 K 3 K 6 


K 4.A 2 7- K 3 JA~A 


(82) 


where 


sin (f - L r ) = 


cos Q 5 




tan (K^ Kg + Kg Qg) 


- VS 


(A - q|) 


Vp = JA (K 2 - Q 2 ) - 2 Kg K 4 K 6 Qg - A (k| + K !) 


cos f = 


A - k ^4 


Q 4 /k 2 - A K 5 


Let 


°2 = 


L 4 (K 1 ~ K 2 } (K 4 - kQ 4 ) K 6 (K 4 K 5 + K 5 ^4 ~ 2k ^ S 


K 4 " uQ„ (k 2 - A k 5 ) " K 4 uQ„ (k 2 - 5 k R ) 


4 5 




(83) 
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Then: 


Q 2 = vG 2 + 


Q 3 sec Q 5 

id 


Vi - A sin (f - Lg) 


cos Qg sec QgGg - 2 ” K 6 sin " L 6^ 


«i* s ^ - «4«i * -" 2 * - v 


(84) 


where 


sin (f - Lg) = 


K 4 u ®4 cos L 6 “ ( K 4 " kC V sin L 6 


®4 -A 2 - K 4 K 5 


One notes on the inspections of G^, and sin (f - Lg) that these are not functions of 

L_, L„, or K . SothatF/m |A (x) ] and H do not contain these variables. This implies 
7*3 ^ • l J o 


that K = 0, K = 0, and L = 0 as desired. 
7 o u 


Formally we may write 


. , ,2 _.2 „2 

a(x) = Qj + Q 4 


tt-2 t^2 

_,2 2 ^7 ' 6 . 2 « \ 

G 2 K 4 + 2 Sin (f “ L 6 


(85) 


We now proceed to write the modified Delaunay variables. 
Let: 


,3/2 


K l = 

_ a 

1 

5 

K 

i 

k 

K 2 = 

a 2 

A* 

K, 

k 

w 

Ol 

II 

a 5 


K , 

-a* 

11 

«7 




3 3 

6 = a 4 
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Then we will have p^ 9 p^ 9 and as before with /3g anci ^4 similar to a ^ and a ^ before. 
The Delaunay -like variables will be j3j, a^, a j8g, and j3^. The relationship for w is 
given through 0^. The new generating function becomes: 


a 


3/2 

T~ “l t + 


.3/2 


« 2 Q7 


+ (a x + or 2 ) 


uQ 4 a 5 /2 -1 k ~ a 5 Q 4 

— + cos 


/i 2~ 
V k "Ve 


- Q l u ' V + a 3 


Qg - sin 


-1 Q 3 tanQ 5 -1 

+ cos 


a 4 a 6 +a 3% 


f a 6- Q 3 7*6 - 7 / 


7" a 3 


+ a, 


cos 


-1 


4 - kQ 4 


Q/ 


, 2 2 

k - a 5“6 


- sin 


-1 


a 6 sin Q s 


2 „2 

a 6' Q 3 


+ COS 


-1 


V6 +a 4 Q 3 


j4-4j4~~ d 


+ sm 


-1 


a l% +a 6 a 3 a 4 


2 2 


a 6/7‘ a 3/7“ a 4 


Our a's are now defined as 
k 


a* = 


a 


372 H o 

5 

k 


372 


(87) 


r (*2 Q 3 sec2 Q 5 tanQ 5 ' a 3 Q 3 sec2 Q 5 " p 2 v) 


(88) 


ffr = 


2k a 6 2 


2 2 „2 2~ 

* 6 = v + Q 3 sec Q 5 

2 ^ = (XgV - C^Q 3 sec 2 Q 5 + of 3 tanQ 5 ) 2 + (p 2 - k 3 QgtanQ 5 ) 2 + 
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where the equations for a g and a g define u and v of the generating function. 
Next we record the (3's: 


/3i = - 


4^ , U % a Y 2 -1 k ' a 5 Q 4 

-4 — t- + cos 


y* 2 - a 5 a 


09 = 


.3/2 J/2 

Qrj + 


k ^ k t + l3 l 


Q 3 tanQ 5 -1 a 4 a 6 + a 3 Q 3 

/3„ = Q„ - sin — + cos 


/ 2 ^2 

V Q 3 


y 2 6 - Q2 3 yp 


0 4 = COS 


-1 a 6 " kQ 4 


Q 4V k 2 " a 5 a 6 


a sinQ _! 

- sm — ■ — + cos 

/ 2 „2 


a 3 a 6 +a 4 Q 3 


/ 2 „2 /2 2 

V Q 3 VV a 4 


3a 


1/2 


0* = " 


2k 


1 1 2 

( a 2 Q 7 \ ( a l +a 2 )a 5 

IV --T“ J + “‘to! 


i q 4 - 


2 /2 , x - 

a 6 (a 6 ~ ^ 

2Q 4 (k 2 - a 5 a 2 6 ) 


Q 1 a 4 a 6 (a 6' kQ 4 ) 
2u 2uQ 4 (k 2 - a 5 a g ) 
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a 3 / 2 (orj + a 2 ) a Q {a\ - kQ 4 ) Qj <* 6 Qg 


kuQ 4 (k 2 - a 5 a g) 


uQ 


Q3 (g 6 a 3 +g 4 Q 3 )tanQ 5 + Q 4 (a 5 g 6 +Q 5 kQ 4~ 2k2) 
v ( g 6 " ^ u Q 4 ^ ' a 5 g 6^ 


Q 3 a/p 
g 6 ^ a 6 " 


fir, = sin 


a 7 Q 3 +a 6 g 3 g 4 


2 /2 : 

a 6v 7 " a 3 >7- 


where 


P - « 7 (ag - Q 3 ) - 2a 3 a 4®6% " g 6 ^ a 4 + a 3^ 


We proceed as before to determine the \ T s needed for A(x): 


3ua 


1/2 


Ql - 2ufi 5 + k 


a l t - 


a 2 Q 7 \ 2(n 1+ a 2 )a^ /2 


Q„ 


g 6 ^ a 6 ~ k< °V 

2Q 4 (k 2 - a 5 a 2 ) 


a 4 a 6 (a 6 ‘ kQ 4 ) 
„ ,, 2 2 X 
Q 4 (k -v 6 ) 


9a " 


Q /« tanQ (a 6 n -i-^Q ) 
vG 2 + "T - 

« 6 V 




v V'P 
(«! - Qo) 


G = . ?§ _ g 5 /2(a l +a 2 ) (g 6 :_ kQ 4 } + _5l + g 4 (g 5 g 6 + a 5 kQ 4 ' ^ 
2 a 6 kuQ 4 (k 2 - ff 5 a 2 ) uQ 2 a 6 u Q 4 (k 2 - a 5 a 2 ) 


VCOSQg 

cos Q 5 Xg = Q 3 sec Q 5 G 2 - — g 


a 6 tan Q 5 (a 6 <* 3 + a 4 Q 3 ) - v VP 


2 f ~2 

a 6- Q 3 


sin /3 P 


a 7 Q 3 +a 3 a 4 a 6 


/2 2 /2 2 

a 6*y a l ' a 3 a/ * 7 " a 4 


COB - 1 - tt6 ' kq » 

%,/ k2 -vl 


sin (f - 0 4 ) = 


( a 3 a Q + ® 4 Q 3 ) a Q sin Q 5 - v cos Q 5 Jv 


, 2 /2 2 ~ 
( a 6 - Q 3 V a 7- a 4 


so that 


and 


cos Q 5 = Q 3 sec Q 5 G 2 - *7 sin (f -0 4 ) 

a 6 

Qo SeC Qjr /~~n n 

Q 2 = vG 2 + 2 / a l " a 4 sin(f - jS 4 ) 

Oin V 


q| + cos 2 Q 5 X 2 = a 2 G 2 + a \ 01 4 - sin 2 (f -0 4 ) 


(91) 


sin (f -0 4 ) = 


« 6 u Q 4 c °s/ 3 4 - (ffg -kQ 4 ) sinj3 4 


Q 4 a/ k2_ a 5 a 6 
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o o o 

Again we note that Qj and Q| + cos Q5 X 3 do not contain a 3 , P7, and P3. Therefore, 
= 0, = 0, and a- = 0, as desired. 


-3 . v 7 3 } o*i\a <~*ry 

We formally write: 


A 2 (X) = Qj + Q 2 


G> a* sin 2 (f -13 4 ) 


( 93 ) 


The same Q’s are used in this problem and the Q is the same as previously cited. 
Below are listed the a's: 


a 3 = 0 


Clrj = 0 


a 


FQ| (vQ 2 + Q 3 X 3 ) 


6 mA(x)a 6 

«5= - 2 mW (vQ 2 + Q 3 X 3 + uQ l ) 
lr TT 3<Xn F 

“ 2 = ^37^ ^ X) + O^mapo- (vQ2 + Q 3 X 3 + uQl) 
5 

.2 

6 


(94) 


f q| r o 

a 4 = m A( x) a2 L 2Q 3 *3 *2 ^ an ^5 " ” p 2 h ~ vX 3 sin Q5 cos< 35 


“o^ ( vQ 2 + Q 3 X 3 ) ] 


. 3 oi^ F 

“1 = a 5 m A( X") (vQ 2 + Q 3 ^ + uQ 1 )+ m ^3/2 |"l 

5 


[l^" (Q 2 v+Q 3 X 3 )+ p l] 


-uQ 4 ^ -2Q 3 XgtanQgH- P 2 + x 3 [ x 3 °osQ 5 (-uQ 4 cos Q 5 + v sinQ 5 ) + p g 


- 3/2° 2 ***> 

n / m 


a 5 m 


We have now presented the information to proceed as before. 

We develop the Poincare variables in two steps. First, we consider the transformations 
of section II and note that we need the following transformations: 
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K*3/ 2 

K 1 - - -i— K 1 

k *3/2 

*2 = -^ Kijj (95) 

Kg = -K| + Kg - K| 


K. = Kf i = 3, 4, 5, 7 . 

The K's with the asterisks are the new constants. Our S function then becomes: 

K 1/2 


S = + 


K*3/ 2 
5 


q 7 \ 

- KJ t + J + ( KJ + K|) 


uQ 4 Kg J 


■ + cos 


-1 k " K 5 Q 4 
'k“ -K^g 


1 Kf-kQ 4 j K|sinQ 5 

-cos -- + sin 


“W* - K 4 ** 


y K i 2 - <y . 


- Qj u - Qg v + K* 


Q 6 - sin' 1 


Qg tan Qg 

y K i 2 -«i j 


(96) 


( _i Kf-kQ 4 _j K*sinQ,- 

+ Kg cos % 4 - - sin 1 4 5 


Q, 




KJ 2 K| 


y K i 2 - ^ 


+ sin" 1 


K 7 S + K 3 K 4 < K 6 ~ K ! - K 2> 


+ (Kg - K| -Kjjj) cos 


K| y K * 2 - K* 2 y K * 2 - (Kg - K* - Kg) 2 
-1 K 3 K 4 + Q 3 (Kj*-*q-Kg) 

y K s -^y^-n-Ki-^ 12 


+ cos 


-1 KSffig-iq-R^RgQ., 

y-i 2 - ® 2 ,/v - k s 2 
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Our K* T s are now defined: 


K i = ^ H o 


*5 " K p/ 2 rX 7 

K 3 = P 3 

K| 2 = v 2 + Q 2 sec 2 Q 5 


(97) 


K* = 2k_ . fi. _ u 2 

5 Q 4 O 2 

K 6 ^3/2 ^ K 1 “ K 2^ + KT ^ ^2^3 sec °5 tan °5 " p 3 °3 sec2 °5 “ p 2 

5 4 

K^ 2 = ^ X 3 v ” *2 ^3 sec2 % + PstanQg) 2 + ^ p 2 “ ^3 Q3 tan Qg) 2 + Kg 2 . 


We next record the L’s: 

,3/2 


L i 


K* 


+ cos 


t - 


u Q 4 K 5 


1/2 


_i K|sinQ 5 
+ sin — - - 


k i 2 - 4 


(98) 


u Q 4 

k+K| 


- cos 


-1 


K|K| + (Kg-Kf-K|)Q 3 


y K l 2 " <4 y®5 2 -(K$-Kf- K *Y 


H- 


K| 3/ 2 q 7 


L l 


K 5 


3/2 


L§ = Q e 


. -1 Q 3 tan Qg _! K| (Kg - KJ - K|) + Kg Qg 


sm 


+ cos 


y k s 2 - ^ /k } 2 - oiyK*, 2 - ki j 


l i= 


COS 


.1 K l - k «4 


. -1 KjstoQj .1 K3Kj*(K5-K|.K|)Q 3 

sin • — — + cos 




y k 1 2 - 4 yKi 2 -^>*, 2 -n-Ki-Kj) s 
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= sin" 1 


K* 2 Q 3 + K| Kg (Kg - K| - Kg) 

K| (Kg-K| -Kg) 2 Kg 2 - Kg 2 j 


3K*^/ ^ q 

L 5 = ZE H -T- - KJ t) + (Kf + K|) 


K| 2 - 2kQ 4 -K|K* 5 1/2 Q 4 
" 2K| 1/2 uQ 4 (k+K|K* 5 1/2 ) 

L| = (KJ + Kg) 


uQ 4 K| 1 / 2 Q 4 
2kK| i/ 2 + 2uk 


KgK|(K| 2 -kQ 4 ) Q x 
2 uQ 4 (k 2 -K| 2 Kg) +_2Tr 


k + KIK^/ 2 (k-Kg Q 4 ) Q 2 tanQ 5 


kuQ. 


u Q 4 (k+ K| K*^) v (K| 2 - Q 2 ) 


Ql K I Q2 K 1 KgK|Q 3 tanQ 5 


uQT 


-I 1 U 

* v (K* 2 - Q 2 ) 


K 6 


k (k _ Kg Q 4 ) x Q 3 tanQ 5 
uQ 4 (k 2 -K| 2 K* 5 ) + ^ ’ v(K* 2 -Q 2 ) 


J*t [ 


Kg 2 - (Kg - Kg - K|) 2 - Kg 2 ] - 2K|KgQ 3 (Kg-Kg -KJ) - K^Q 2 G> 3 
K| (K| 2 - Q 2 ) 

L * _ Q 1 (k2 + K 4 2 K t + kK 4 K *i / V Kg K* 5 3/ 2 (K| 2 + k Q 4 ) u Q, 

4 (Q| K* 5 3 / 2 - K| k) (k+KIKI 1 / 2 )" ’ (K| 3/ 2 Q 2 - k K|) (k 2 - K^Kg)" 


, P 1 Q 4 

(K* 5 3/ 2 -K|k) 


l-, r , r ,l/2 , k < k - K W -2^^% 

4 5 k+K^Kk 1 ' 2 _ K| (KJ Z - Q^) 


Kg v Q 3 tan Q 5 Qg (K| 2 + Q 2 tan 2 Q g ) v 
K| (K| 2 - ©|T 


K| (K| 2 - Q 2 ) 


Q 3 Jk I 2 [k * 2 - (Kg - KJ - Kg) 2 - K| 2 ] - 2 KJ Kg Q 3 (Kg - Kf - Kg) - K* 2 Q 2 

9 _9 


K| (K* 2 - Q 2 ) 


53 



L § 


3K* 5 


1/2 


2 k 


/ K+ Qt. . _ K|K|uQ 4 K| 1/2 ( K | 2 + kQ 4 ) 

^ <r 1 2 (k Z - K* 2 Kg) (K* 3/ 2 Q 2 - K*k) 


*4 


uQj 

2 K* 5 L/ ' 2 (K* 3/ 2 Q 4 - K| k) 


v 2 „2 . K f k 


-K|Q3-Kf + 


(k+ K| K^ 1 / 1 ) ; 


hQl / o (K| 2 -kQ 4 )k \ 

+ 1 /» 0-7^9 n ( - u Q d + KJ Q, + k t7 o I 

2K* 5 1/a (K* 3/2 Q^-K|k) \ 4 5 4 Q 4 (k + K|K*f rj ) / 

9 2 

One notes that the factor k -K| Kg still remains in Lg, and LS. To remove these 
we need two transformations. One we rename Lg and Kg. The other will affect Eqs. (46) of 
section n. This will not be done directly because of the needed constants. These trans- 
formations , which constitute the second step, are 

S = Lrj£ Lg + identity 


_ •£>-* __ y T 

Trqr ~ ^6 


8 S 


3 L, 


_= - H = L§ 


* -/ k - K 4 S V2 


sin K 6 ^6 


+ / k - K i K 5 


1/2 

5 COSK 6^4 




as 

9K' 


as 

aic^ 


= - Li = 


i = 1, 2, 3, 7 
[+ sin Kg/3 g - cos K£/3 4 ] 


2 Jk - K 4 K g 


1/2 


+ K 4 [sin Kg /3g - cos Kg 0 4 ] 


L 5 = 


4 k; 


1/2 


> - K 4^ /2 


- 0 P 


as 

9K 6 


= - L 


I 1/2 

g = " i k - K 4 Kg (cos Kg /3g + sin Kg /3 4 ) 




(99) 


as 


-=~ a 4 =1^- 


1/2 

K 4 Kg cos Kg 
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as 


W = -0!6= ' 

W = _Q! 5 = " K 5 


k-K'K’ 
4 5 


1/2 


sin K’ 


The other transformations result in identities. Our transformations are: 


K i =«i 


h =?i 

K 5=“5 


i = 1, 2, 3, 7 


. 2 2 
_ k " a 4~ a 6 
K 4 --^- 


H = - J k " 


* * 1/2 

cos Lg 


K 6 " "“ 4^6 + a 6@4 


, * 1/2 
a c = - | k-K 4 K* 


sin L£ 


L 1 


1/2 

( a 6@6 +a 4 < 3 4 > 


2(«4 + a |) 


L*= - 


( k -^-« 2 ) ( a 6 / 3 6 + « 4 / 3 4 ) 

- rt n + PC 

4 "5 


1 

L£ = tan™ 1 * -5- 

6 a 4 


Pa = 


_ 2L| /k -K|Kg 1//2 cos Lg Kjji sin Ljjj 


K* c 


172 ^ 


k-K|K* 




+2 L*/k-K*K* sin L* K* cos L* 
^6 ' T7-+1/ 2 


K r 


/ k - K| K* 1/ 2 


K | 

^5 _ L 5 " 2 K* L 4 


( 100 ) 
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Next we will rewrite our equations in the new terms. 


k H 


*2 = ~~372 


a 3 = p 3 

2 k 2 v 2 + Q|sec 2 Q 5 

= o^- u g- 


a 2 = (Kg v-XgQgSe^Qg + PgtanQg) 2 + (p 2 - K 3 QgtanQ 5 ) 2 + a 2 

Of 2 Qg+ Pg (>^ Qg sec 2 Q 5 tan Q 5 - Pg Q g sec 2 Q 5 - P 2 v) 

Sm ^ 7 Ja\ -a\ Ja\ (v 2 7 Q 2 sec T Q 5 ) - o| sec 2 QgtanQg - P g Qg sec 2 Q 5 - p 2 v) 


v 2 + Q 2 sec 2 # -kQ^ 


k+ aV 2 / v 2 + Qo sec Z Q 5 


-=J k - 2 yC 2 + Q 3 sec 2 ^ cos f 


sinfl 


cosH 


4 2 + Qg sec 2 Q 5 J a 2 q - a\ cosJ3j 
a 7 /v 2 + Qgtan 2 Q 5 


(v 2 + Q 2 sec 2 Q 5 )a 3 + Qg Q 2 sec 2 Q 5 tanQ 5 - P g Qg s^Qg-^v)^ 
it 2 + Q 2 tan 2 Qg^/a 2 (v 2 + Q 2 sec 2 Q 5 ) - Q3 sec 2 Qgtan Q- - P g Qg sec 2 Q 5 - P 2 vY 
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4 ' 2 «7 


+ 01 + 


/3 3 = Q g - sin 


_! QgtanQg 


4 2 + Q^tan 2 Q 5 


(k-ag-of ^) 2 2 2 2 

. -ij a 7' K 6 K 4 

+ sin N 


“ 7 y v2+Q 3 tan2Q 5 


2 2 2 
K 6 K 4 = ^2 ^3 sec Qg tanQg - a g Qg sec Q 5 - P 2 v 


, 2 2 
k-a 4 -a 6 


[(« 6 0 4 -‘W " ( a i 


'6 a 4 2 2 

a 4 + 01 6 


~ a 4 a 6 


where Kg and two forms of L| are given in Eqs. (97) and (98). 

„ 3 4 /2 «T . „ . 'iS 


< a 2 — “ ff l “ 2a7~ + 


2[(k-a 2 -a 2 ) 2 -a 5 Q 2 ] 


/(k-ag-a|) o 2 \ 

Q 2 v ( «5 +Q 3tan Q 5 J 


P 2 tan Q 5 - a> 3 v Q 3 tan Q 5 


yj [ (k " a 6 - a 4> -“5^][°7 ' ( V4 • Ve ■ a l • a 2 )2 ] - [ a 3 O' -®6 - a 4 ) +Q 3“5 1/2 (®6' 3 4 -°4 (J 6 '“l -»2>] 
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ft* = 


3 “5 /2 , % 


2 k 


(a 2~ -"l t)+ 2^7l Q 


a 2 5 Q 2 -( k - a 2 6 -^) 2 

a 5 uQ 4 


s 


2 9 2 

(k-«6-<) 

a 5 v 


“3 (k ’ “6 ‘ “4^3 ^ Q 5 , ( “6^4 * a <P 6> (k ' a l ‘ a 4> 

2 

Q^Ogtan Qg 

r 2 2 2 2 "I a 1 / 2 

V f k -V a 4 } " a 5 Q 3 J 5 

uQ 4 ' 

v [( k ‘“6'“4 ) '“5^] 


+ («!+ « 2 ) 


2 (o-g Q 2 - k (k - a 2 - a 2 ) (k - a 2 -,a 2 ) a^Q 2 tan Q g 


kuQ 4 a 5 


L 2 2/ „2l 

• (k-a 6 -a 4 ) " “5 Q 3 j 


+ Vs /2 J 

(k-a 2 6 -* 2 4 ) 2 -a 5 Q 2 


”7-<¥4- 

Ve ~ “i " °a> ] - 

[ a 3 < k - *6 ' a 4> + % a 5 2 <^4 a 6- 


<k-« 2 -* 2 ) 2 -« s Q 2 ] 



One obtains the following time derivatives: 

®5 = ‘ H^y (uX l +vX 2 +Q 3 V 

• _ 3k ^ •. , k(r F 

“2 “ m 2 <*> 


« 3 = 0 


a 7 = 0 


FQ, 



X 3 1^7 co® 2 V P 3 (2^Q 3 tanQ 5 - p 3 -vX 3 sinQ 5 cosQ 5 )J - P 2 P 3 ^ 


■]- 


cos/3 7 mA(x)ya 7 - « 3 „/a 2 (v 2 + Q|sec 2 Q & ) - |( >^Q 3 tanQg - p g Q 3 ) sec 2 Q 5 - p 2 vj' 


* e 3 


stn<3 7 V^-°i{°7 (v X 2 +Q 3 V P 3‘ V ^!i nQ 5 COsQ 5 ) •_ P 2_ X ?]}] 

{°7 (v 2 + Q§ sec 2 Q s ) - [(^Q 3 tan Q 5 - Pg) Q 3 sec 2 Q 5 - p 2 v] ]• 
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* F Q 4-v/“7 ~ a 3 cos ^7 v ^3 v cos2 Q 5 ' X 2 Q : 


H = 


* 3 ) 


- tan H tanP _P 


mA (X)a 7 (v^ + tan 2 Q 5 )^ 2 ^ / / v2 + Q 2 sec 2 0 & cos H 


7 7 


Tf< 

1 

1 f F «4 

/ v 2 + Q 2 sec 2 Q 5 -kQ 4 

cosH sinH 

1 

CD 

l».y° _p* Q | „ c 2 q 5 jv 2 ^ 2 ,»» 2 « 5 r i<x) 

V Q * 

sinH -cosH 




2 „2 2 _ 

+ Qg sec Qg u 


-sinH cosH 
cos H sin H 


'-2 F Q 4 (X 2 v + Xg Q g ) 

‘ — rriA - (x) 


x g cos Qg 

(v x 2 + Q 3 k 3 ) sin Q 5 


^/v 2 + Qg sec 2 Q E 


+ A /v 2 + Q 2 sec 2 Qg u Hj 


cos H sin H 
sin H -cos H 


V 2 + Q 2 sec 2 Q 5 - kQ 4 . F [u Q 2 (v\ 2 + Q 3 Xg) + \j (v 2 + Q 2 sec 2 Q g ) j 

Q, H - 

mA 


(x)^/v 2 + Q 2 sec 2 Q 5 


- sin H cos H 

cos H sin H 


v cos Q c 


/ 2 _2 2 _ . _ 
v + Qg sec Qg sin Q g 


-( 


FQ 4 (v>^ 2 + Qg X g sin ^ Qg) \y / v 2 + Q3 sec2 Qg ( k + + ^3 Sec2 Q5 ) 2a 5^ 2 + ( y2 + Q3 ta ° 2 Q5) “5] 

mA (X) 2a^ /2 ^/v 2 + Q 2 sec 2 Q 5 ( k + ff* / 2 ^/v 2 + Q 2 sec 2 Q 5 )(v 2 + Q 2 tan 2 Q & ) 


Define 


“5V^ 2 + Q 3 sec2 Q 5 


4 4 ^ (^^ / Vv 2 + Q 3 2 sec 2 Q 5 ) 


a 5 /2 V^ 2 + Q 3 ten2 Q 5 

~ 2 2 

k - a R - a . 

b 4 


sin 1 = 
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H 


CO 

<0 

v. 

f_,. 1/2 . . . . 

FX l a 5 | ( a 6 a 6 +Q W u 

fq 4 

(u 2 -a 5 )^i + (v ^2 + u 

(a 5 Q 4 +k-ag-a 4 ) 

mA(x) + (2 k -4 - 04 ) 

v. y 

1 mA(x)ag^k 


FQ|- cosi ^ " sin ^ Q5 ^3 cosQ^) tanQg cos i 3a^\ 

mA(x) sin iVv 2 + Q 2 tan 2 Qg 


5 2k 


•\/ a 2 - “3 i2 ( a 6 * 6 + « 4 * 4 )cos^ 7 + (k -a 2 -a 2 ) 


• /® 5 FQ^(vX«+Qo^oSin^Qi-) 

sin^ 7 /3 7 + cos^ 7 f— + - 3TJ 


"5 mA( x )( v ii + Q 2 tan 2 Q t -)/J 


v{k -a| -a 2 ) (a 7 sin 2 |3 7 + a 2 cos 2 |3 7 ) - » 5 a 2 Q‘ 


3aJ /2 a. Q 7 
5 5 


h ~ 2¥ 

^3 = 0 


= ®i + *2 


( K g- a i- a2 )FQj( v Q2 + Q 3 V 

mA(X) (v 2 + QgSec Z Q 5 ) 


FQ 2 [(2Q2 Qgtan Q 5 - a 3 - v sin Q 5 cos Q 5 ) X3 - ^Sa] 


mA(X )>/ v 2 + Qo sec 2 Q r: 


a l " “ 2a7 


3 «1 k r 


F 

mA (x) 


_K 1 ( ^( vQ 2 + Q 3 X 3 > + P 1 ) 

_X 2 (-uQ 4 X 2 -2Q 3 XgtanQ 5 + p 2 ) + X3 (- XgCosQ 5 (vsinQ 5 -uQ 4 cosQ 5 ) - pg) 

+ —^2 A(x)} 

m J 

if 

The time derivatives of L 4 and /3 c may be obtained in a similar manner. They will not be 
presented here because of their length. 



I 


The following equations give information to define A(x): 


+ (a i + “ 2 ) j -rm ~ 


u 2 Q 4 


a \ n % 


(k-a|-a|) 2 

2 kQ^ - (k - afg - a 4 ) 

«J /2 Q 4 (2k-«|-a|) ' 


K*K| 




/ (k -a? -a 

(- 1 :- 


Q 4 ( a 6 + a 4^ (2 k- Q!g - « 4 ) 

where 

K|Kg 

(« 1^5 

is well defined except at ag =a^ = 0. 


K|K^ 


*1 + a 2 


(a 6 +a 4 ) a g ' 


2/3 5 l; 

of ’ Q$af*(k- a %- a %) 


( Qfj + a 2 ) v u [(k - a| - a 4 )“ + a 5 k Q 4 J 3 a ^ 
ag 2 Q 4 k (k - ag -aj) (2 k -aj| - a*) ' “^k 

K|K|vu|i-a 2 -a 2 ) 2 + a 5 kQ 4 ] <\ 


(«2 — “ a l fc ) 


N v a g 
„,2” “2X 


„ 2 2\2 /-> / 2 277,7, 2 ZT ,, 2 272 r„ 2 Z72 73 

(k-a 2 -a 4 ) Q 4 (a 6 + a 4 )(2k-a 6 -a 4 ) (k-a 6 -a 4 ) (k -ag -a 4 )-a 5 Q 3 


Q3 ta-nQ^ag^ 


7, 2 ZTTh 2 272 77211 T 

(k-a 6 -a 4 ) (k-a 6 -a 4 ) -a 5 Q 3 J a g 


a 3 (k-a| -a|) 


+ Q 3 (K 6- a l- a 2 


I 



where 


N = 


(k-a 2 -a 2 ) j® 7 " ( K 6 -a i" a 2 ^ J"j^ Q! 3 ^ k " a 6 “ Q! 4 ^ + % a 5 /, 2 ^ K 6 “ Q! l 

2 2 r 2 ^5 L5 Jag Q 4 - (k - a„ - a ,) 

sec 2 Qg tanQg i * - - 4 - l 5 ^ _6 4 

3 5 HOT Q 4 < (k-a 2 -a 2 ) 


,V2 


Q 4 k 


Q 7 

(a 2 T“ 


(a 1 + a 2 )u[(k-a 2 -a 2 ) 2 + agkQ 4 ] 

" + « 5 / 2 ( k -« e Ta 4 ) Q 4 k ( 2 k - “ I " « 4 > 

Kg u [(k - a| - a 2 ) 2 + a 5 k Q 4 ] a/n a|/J Q 3 

(k - a 2 -a^) Q 4 (a 2 + a|) (2k - a 2 -a 2 ) (k -a 2 -a 2 ) 2 ^ [(k -a| -a 2 ) 2 -a 5 Q^] 

a 3 Q 3 sec 2 Qg v a 5 (Kg - aj - a 2 ) [ (k - a 2 - a 2 ) 2 + a g Q 2 tan 2 Qg ] v a 4 / 2 

(k -ag -a 4 ) [(k -a| -« 4 ) 2 -a 5 Q^j 


Pi = 


(k-a 2 -a 2 ) 2 -a 5 kQ 4 

5 5 ?- *-1 2j3g + 


“5«4 


{ 


(k-a 2 - a 2 ) 


1/2 


3a: 7 “ Q 7 

L I -- 4 — <“2 4 ““I*) 


a 


3/2 


5 


- (a 4 + a 2 ) u 


a 


, 1/2 


/, 2 2,° 

(k - a 4 - a g ) 

. 3/2„2 ~Z 2 2( 

ka 5 7 Q 4 ( 2 k - a 4 - a g ) 


[ 2^5 L* [a 2 Q 2 - (k - a 2 . a 2 ) 2 ] 3 a* /2 Q 

x 0 =< - — (a 0 - a+ t) 

d 1^2 _ 9. 3/2, 2 2. .2 2 cr 1 


* * 9 9 2 

| K 6 K 4 u(k- a 2 -a 2 ) 
a 5 (^ (a 2 +a 2 )(2k-a 2 - 


Q 4 k 


~2 3/2,, 2 2, 

Q 4 Q 4 « 5 (k ~ a 6 -a 4 ) 

V 

(a x + a 2 ) u [(k - ag - a 2 ) + agk Q 4 ] v/N a j/ 2 Q 3 

■( /ft n « ft ft _ . r 


a 5 /2Q 4 k(k ' u 6 


4 7 + “5 K ^4J v " “5 ^3 

2 2 , 2 2 ,,, 2 2,2 f„ 2 2,2 21 

-a 6 -a 4 )( 2 k-a 6 -a 4 ) (k - a g - a/ [(k-a g -a 4 ) -a 5 Q 3 J 


K 6< u[(k - a a 6 . a ^, a5t Q 4 ]] , + -^°5 /2 

0 9 9 9 2 9 9 I 3 Sr 9 o 


Q 4 (“g + a 4 ) (k ' a 6 " a 4 )2(2k _a 6 " a 4 


[(k - a 3 - a 2 ) 2 -a 5 Q 
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tan Q 5 [ a 3 (k - a l - a 4 } + 4 /2q 3 (K 6 “ a l - a 2 ] ] va z 


[< k “ “ «4 )2 " "5 Q 3 ] 

If we proceed as before we write: 

2(k -4-4^5 L 4 [ «| q| - (k - - «^) 2 ] 


A 2 (X ) = q2 + Q 2 




<4 


[( k - «e ■ a 4 } +a 5 kQ 4l 

(<*! + a 2 ) (k - a 2 - a 2 ) - ka^ 2 

K 4 K 6 

(<V*4>- 

L J 

- 


(108) 


a 5 Q 4 k(2 k - a 2 - a 2 ) (k - a 2 - a 2 ) 


( a 2 ~T~ - «!*> < k * a l ' a \) 


kQ 2 


„ 1n (f x„ -1 a 6- 1/2 V a 7 " (K 6 " a l ■ a 2 > 
sm (f- tan — )o' 5 / „ s 


i 2 2 

k -«6‘ a 4 


It can be seen, from the equations for /3j, /3g, and Q 4 that a(x) does not contain 0g, /3 7 , or ttg. 
Again we write: 

3/2 


H = — — a + F ^W- 

XX . w, T r- 1, 


-372 (^2 " ^1) ~ * 


and we obtain 


a 3 = 0 


= 0 


^3 = 0 
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APPENDIX 


The following tables of equations were found useful in developing the equations presented 
in the text and in developing computer programs for numerically evaluating the ordinary dif- 
ferential equations. 

Only equations for the initial constants (K, L) of section II (K-D-P- variables) and some 
tables for the Poincare variables are presented. Tables for the other variables may be de- 
duced from these tables or derived in a manner indicated by the tables. 

In short, these tables are presented only to aid the interested reader in understanding 
the results and extending them for his special purposes. 



TABLE FOR K-D-P VARIABLES 


Z = -Q 0 


sin i cos (K^ — f) 
1— sin 2 i sin 2 (K^— f) 


+ Qg sin" 1 


+ Qg sin 


cos i sin (K^-f) 


L . 2 . . 2, v 

1-sin l sin (K^-f) 


sin i = 


/ K 2 k 2 
K 4 K 3 
K. 


31 _ 

3K„ 


3i 


cos 1 


ft 2. 
3 l 


K. sin i ’ 3K„ K. sin i ’ KJK, 


a 2. 

3 1 _ _ COS 1 
3K 2 K 2 sin 3 4 i 


A 

3K 2 


cos i (2 — cos 

2 .3“ 

K 4 sin i 


cos f = 


^-X«4 


Q t/ 


k 2 -K,K^ 
5 4 


^sin i sin (K^-f)j 


1 

~ 2 . 3 . 

K, sm i 
4 

i) 


8 , _ ^ K 5 K 4— K ^4 a , -K 4 (K^-tQ 4 ) 

9K 4 TQ 4 (k 2 -K 5 K^) ’ 8K 5 2uQ 4 (k 2 -K 5 K 2 ) 

3f_ = fl 9 2 f _ - 2K 4 K 5 k ^ k - K 5 Q 4 ) K 4 (2k 2 -K 5 K 2 -K 5 kQ 4 ) 

8Q 4 uQ 2 ’ 3k| uQ 4 (k a -K 5 K 2 ) 2 v^Q^^-KgK 2 ) 

3 2f - K 4 (K 2 -kQ 4 ) 

3k| 2u Q 4 (k 2 -K 2 K 5 ) 2 4u 3 Q 4 (k 2 -K 5 K 2 ) 

a 2 f k 2 (K 2 - kQ 4 ) 2k 2 - K 5 K 2 - K 5 kQ 4 

4 5 uQ 4 (k -K 4 K 5 ) 2u Q 4 (k - K^) 
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„ K. cos i cos (K 7 - f) cos i sin (K„ - f) 

- Qo r - — o o — v * 7 -0 + Q, 


9i 


[' 


2 2 
1 - sin i sin 


- Op 


9 2 Z - Q 
9i 2 " Q2 


( r 7 - f)j r y / 1 _ gin 2 . gin 2 ^ 

sin i sin (K^ - f) cos (K^ - f) 

£ 1 - sin 2 i sin 2 (K^ - f)j 

sin i cos (K^ - f) [ cos 2 (K^ - f) - 2 cos 2 i sin 2 (K^ - f) j 

f, . 2 . . 2 ,71 5/2 

1 - sin i sm (K 7 - f) 


- Qs 


sin (K^ - f) cos 2 (K 7 - f) sin i sin (K ? - f) cos (K^ - i) cos i ^ 1 + sin 2 i sin 2 (K ? - f)J 


j~l - sin 2 i sin 2 (K ? - f) J 


3/2- 


[ 


2 2 
1 - sin i sin 


T 2 2 

0 2_ K. cos i sin (K„ - f) cos i - 2 sin i cos 

8 z _ . r> Z 1 

8i8(K 7 - f) ^2 


»7-«] 

2 (K 7 - £)] 


r, . 2 . . 2 5/2 

1 - sm i sin (Krj - f) 


+ <4; 


cos i cos (K^ - f) 


1 - 


sin 2 i sin 2 (K r 


-f) 


372 ' Q 6 


[ 2 2 2 1 
cos (K ? - f) - sin (K ? - f) cos i I 


T 


2 2 
1 - sin i sin 


«,-!>] 


8iaK„ w 2 


cos i cos (K 7 - f) 


, . 2 . .2 

1 - sin i sin 


(k 7 - f)] 


TJ2 


2 

sin i cos i sin (K^ - f) 

8 (K 7 - 1 ) = + Q 2 U _._2 7 ~2 ~ + Q 


dZ 


sm l cos 


(K.J - £) 


Ti . 2 . . 2 «\1 3/ 2 o y" ~ “ « 

[1 - sm 1 sin (K 7 - f) J /l - sin 2 i sin 2 (K ? - f) 


+ Qc 


COS 1 


j^l - sin 2 i sin 2 (K 7 - f)j 
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a 2 z 


a(K ? -if 


= Q 2 


K 4 sin i cos 2 i cos (K ? - f) [l + 2 sin 2 i sin 2 (K 7 - f)j 

5/2 


[l -sin 2 i sin 2 (K^ - f jj 


-sin i sin 
L 

_ 2 . . . 2 . 


sin i cos i sin (K„ - f) 2 sin^i cos i sin (K„ - f)cos (K„ -f) 

• r~\ i 7 

[*- 


‘ Qs [l - sinh sin 2 - '^ + Q « 


(K ? - f )] 


2 2 
sin i sin 


(K 7 - f)] 


a 2 z 


2 

0 ^ sin i cos i sinOEC^ - f) 

3(K 7 - f » 8K 4 ' ^ [l - sm 2 i 3l n 2 (K 7 -'f)']V2 


8 z sin i cos (K 7 - f) 

^4 = r 2 2 

V 1 - sin 1 sm (K^ - f) 


We note: 


az 


az ai az aK 4 

+ — + 


az 9 (k 7 - f ) 


aq- ai aq. 9K 4 aq. + a(K 7 -f) aq. 


9 2 z _ a 2 z ai ai a 2 z f ai 9K 4 ai aK 4 

9q i 8q k ai 2 aq j aq k + aiaK 4 V 9q j 9q k + ^k ^ 


az a 2 i 


H) ai 9q j 9q k 


a 2 z 

(bi 8 ( K 7-» 

ai 

3(K 7 -f)> 

\+ __3 2 Z | 

( 3(K ? - f) 

8K 4 

3(K 7 - 

f) 

3K 4 X 

8i d(K^T) 

l 8q j ^ + 

Sq k 

8q i / 

1 8K 4 8(K 7 -f) 1 

V 3q j 

8q k 

8q k 


*h) 


a 2 z a(K 7 -f) 8(K 7 -f) 


0(K 7 -f) 2 9q j 


az 9 (K 7 -f) 




+ 8(K ? - f) 8q 8q fc 


where q. and q fc are Q 4 , Kg, K 4 , K g and K 7 . 
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Z(3) =cos 


•1 k - K 5 Q 4 
J k 2 - K^Kg 


8Z(3) = _ K 4 (k - K 5 Q 4 )K 5 / 2 
8K 4 uQ 4 (k 2 - K^Kg) 

BZ(3) ^ VV ^4 ~ ^ 4 

8K 5 2K^ /2 uQ 4 (k 2 - K 2 K 5 ) 

k 1/2 

8Z(3) = 5 

9Q 4 uQ 4 

9 2 z( 3) _ ~ K 5 /2(k - W [ u 2q 4 (k2 + K 5^ + n 2 ^ 2 ; k 4 k 5 )] 

8 2 K 4 " u 3 Q 4 (k 2 - K^Kg) 2 

8 2 z(3) _ k 2 K 4 (k - K 5 Q 4 ) Wfe + **4 - 2k2Q 4) 

9K 5 9K 4 ~ " K 1 / 2 uQ 4 (k 2 - K^Kg) 2 2u 3 Q 3 K* /2 (k 2 -K^Kg) 

8 2 Z(3) = K 4 K 5 /2 

8Q 4 9K 4 u 3 Q 3 

3 2 Z(3) ‘ (I 4 - kC *4> k + ‘V*’ - K 5^> ^ “<*4 - k Q 4 K k2 -31^ %) _ t 

4Kj /2 u 3 Q 4 (k 2 -K 5 K^) + 4K 3/2 uQ 4 (k 2 - K^) 2 


8 2 Z(3) 

9K 5 9Q 4 


2kQ 4 - 


A 72 , 
^5 


2Kc 7 “u Q‘ 
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TABLE OF POINCARE VARIABLES 





+ a 4 ) (2 


k 





4 


9 i 
9a 6 



9i 

9a 5 


7 

8 “3 




, 2 2 , 

(« 3 + « 4 ) 
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9K 7 . 

“6 92 

K 7 

- 2a 3 a 4 

9 2 K 7 

a 3 * 


If 

2 2 ’2 
a 6 + a l 9a 3 

, 2 2-2 ’ 
(“3 + a 4> 

9a 3 9a 4 

' <*3 + 

2-2 

*4> 

9 2 K 7 

2a? 3 a 4 

8 2 K 7 

2a 6 a 7 

8 2 K 7 

2 

a 7 - 

2 

' a 6 

9 *4 

,2 2-2 * 
0*3 + a 4 ) 

9a 2 

~ , 2 2-2 
(a g + a 7 ) 

> 

9Q!g 9ff 7 

, 2 
(«6 + 

2,2 

® 7 ) 

0 2 k 7 

“ 2a 6 a 7 







, 2 2-2 
(a 6 + a 7 ) 







<*-»«■ 

2-2 
- a„) 

“ “5^4 





cos f 


8f 

da 


6 


a 5 Q 4 J4 + 4 (2k ‘ “6 ' A ] 

+2ag[k(k - a 5 Q 4 ) + (a 2 + a 2 ) (2k - Og - a 2 )J 

a 5 ^ 2 + “7) ( 2k " a g - “7) 


8f 


+2a 7 |k(k - a 5 Q 4 ) + fag + "7) ( 2k " a g " “7)] 


a 5^ 2 u< ^4^6 + ®7^ ^ 2k “ a 6 " a 7^ 


8f 

to c 


+ (k - <*g - a 2 ) 


■372 


uQ, 


9 2 f 

to 2 


8 2 f 

Ba 6 8a 5 


- 3 (k - ag - a 2 ) (k-ag-a 2 ) 

2 uQ 4 a 5 /2 


/1 2 2 X 2 2 ^2x1 

(k-ffg-a 7 ) -a 5 Q 4 )J 

77273^3 ~ 


2a- “ u Q' 4 


-2a 6 

3/2 “ 
“5 uQ 4 


„ „ 2 2,2 

2ag (k - a g - a 7 ) 

5/2 3„3 
a 5 uQ 4 


9 2 f 

9 7 9a 5 
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uQ, 


O /, 2 2-2 

2a 7 (k - a g - a ? ) 

5/2 3„3 
« k u Qa 


72 



+2k(k - a 5 Q 4 ) 2k(a 2 - a 2 ) + (a 2 + a 2 ) (3a 2 - a*) 

a^ 2 uQ 4 (ag + a 2 ) 2 (2k - a 2 - d 2 ) 2 

_ 2a 2 (k - a 2 - a 2 ) 2 

u 3 Q 4 (ffg + a*) a 5 (2k - a 2 - a 2 ) a g /2 uQ 4 


~ 3/2 3„3 
“5 u Q 4 


sin (K 7 - f) = 


-sin (<i> - Kg) 

/ 2 2 2 ~ 

/cos i + sin i sin ( <f> - 


cos (K 7 - f) = 


cos i cos (<f> - Kg) 

y/ ~2 2 . ^2 

cos i + sin i sin (4> - Kg) 


1 - sin^i sin^ (K 7 - f) 


2 . 

COS 1 

2 2 2 
cos i + sin i sin (cj> - Kg) 


9i ~ " Q 2 


2 2 f 2 2 2 1 

(k - a g - a 7 ) cos (<|> - Kg) jcos i + sin i sin (<4> - Kg) 


1/2 

ffg COS X 


-QgSin (<j> - Kg) + Qg tan i cos (<j> - Kg) sin (4> - Kg) 


9(K„ - f) " Q 2 


22 r 

(k - ffg - a rj) sin i sin (<t> - Kg) 1 1 


2 2 2 
cos i + sin i sin 


1/2 

a g cos i 


+Q 5 sin i cos (4> - Kg) + Qg 


2 2 2 
cos i + sin i sin (<t> - K R ) 


-gg- = -Q 2 sin i cos (4> - K g ) 
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3(K 7 *» f) 


[ 2 2 2 1 

cos i + sin i sin ($ - K fi ) r r 9 9 9 

2 = 2 Qo K 4 cos ($ ~ K g) cos i + 3 sin i sin (4> - Kg) 

, - f) cos i L 1 J 


+ Qg cos i sin (45 - Kg) - 2 Q g sin i cos (4> - Kg) sin i (4> - Kg)] sin i 


a 2 z 


3(K 7 - f)8 K 4 = ' Q 2 


2 2 2 

sin i sin (4> - Kg) (cos i + sin i sin (4> - Kg) 


a 2 z 

S(K 7 - f)9 i 


2 2 2 

cos i + sin i sin (4> - Kg) 

— ' 2 . 

cos 1 


|q 2 K 4 sin (4) - Kg) [3 sin 2 i cos 2 (4> - Kg) - l] 


+ Q c 


+ Qg cos i cos (<t> - Kg) - Qg sin i [l - 2 sin 2 (4> - K g )] j 

- Kg) |cos 2 (4? - Kg) - 2 sin 2 (4> - Kg)j (cos 2 i + sin 2 i sin 2 (4> - Kg) 

1/22. ‘ 
ffg cos 1 

(4^ - Kg) cos 2 (43 - Kg) sin i sin (4> - K g ) cos (4> - Kg) [cos 2 i + 2 sin 2 i sin 2 (4> - Kg)] 

4- On 1 ■ - ' 1 • n 


2 2 

a 2 z (k - ffg - a ) sin i cos (4> 

— = Q “ 

ar 2 


cos i 


2 . 
cos 1 
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TABLE OF NOTATIONS 
Mixed Variables 


sin 0 = -sin (K^ - f) sin i 


y r 2 2 

1 - sin (K^ - f) sin i 

sin Kg cos (K ? - f) - sin (K^ - f) cos i cos Kg 

sin <j) = ■ — -- ■ 

A - sin 2 i sin 2 (K^ - f) 


COS <t> = 


sin Kg cos i sin (K^ - f) + cos Kg cos (K^ - f) 


'l - sin 2 i sin 2 (K ? - f) 


w = K r 


sin i cos (K^ - f) 


/ 2 . 2 

f l - sin i sin (K^ - f) 


sin 0 = 


+ sin (<j> - Kg) sin i 


/ 2 2 2 
'cos i + sin (4> - KJ sin i 

b 


cos 0 = 


COS 1 


/ 2 . . 2 , \ . 2 . 
'cos i + sin (<j> - Kg) sin l 


tan 0 = tan i sin (<f> - Kg) 


75 



= r = 


K 2 ^/cos 2 i + sin 2 i sin 2 ($ - Kg) 


k y4os 2 i + sin 2 i sin 2 (<|> - K g ) jcos i cos K^cos (<t> - K g ) - sin K ? sin (4> - Kg)j- 


v = K 4 sin i cos (<(> - Kg) 

^/k 2 - KgK^ sin cos i cos (4> - Kg) + cos K^ sin (<f> - Kg) 


u = 


K, 


y f 2 2 2 

cos i + sin i sin (4> - K^) 

Poincare Variables 



2 2 v 

“ a 3 ~ a 4 


sin = 


a„ a_ - a . a „ 
3 7 4 6 


/ 2 2 , , 2 27 

/<*3 + ® 4 ) < a 6 + a 7 } 


cos K n 


a 4 a 7 + a 3 a 6 


/ 2 2772 2 , 

y(“3 + “4 } a 6 + a 7 ) 
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sin (K ? -f) 


cos (K 7 -f) = 


a 


1/2 


-(a 


(k-a 2 -< 


°/ 2 <- 

^ + “ 2 4 ) 

(k - a 2 - 

— r 



, 2,2 

*7 “"5 


2 2 

k - a - a 2 

a k (2 -a 3 - <* 4 ) (<* 4 cos4> + a^sin*) 


2 2 
f 6 


a 


Yj2~ ff 3 ~ a 4 ) (<* 4 cos 4> + a 3 sin4>) 


_ 5 
2.2 


k — a 6 “7 2 2. 


/, 2 2 
/2k - a 6 - « 7 


op 1/2 

(k -Qfg -a 7 ) (a ? sin <t> - ciscos 4>) - ( a 3 a 7 - a 4 ff 6 ) ( a 4 C0S <t> + a 


_ a 2 _ a 2 

(k - a J -« 7 ) /(k - - “ 5 (2 T 72 ' 1 " a 3 4 } (a 4 cos 41 + a 3 sin ^ 
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